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^ ■ Abstract: We present a manifestly duality covariant formulation of the composite non- 

QQ ■ BPS and almost-BPS systems of multi-centre black hole solutions in four dimensions. The 

'^ . method of nilpotent orbits is used to define the two systems in terms of first order flow 

• , equations that transform covariantly under the duality group. Subsequently, we rewrite 

f^ \ both systems of equations in terms of real, manifestly duality covariant, linear systems of 

cn I Poisson equations. Somewhat unexpectedly, we find that the two systems are naturally 

described by the same equations involving space dependent abelian isometries that are 
conjugate to T-dualities by similarity transformations. 

X 
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1. Introduction and overview 

The structure of black hole solutions in the supergravity effective description of string/M- 
theory compactifications has long been a useful tool in understanding their microscopic 
realisation in string theory. In the supersymnietric (BPS) case, the supergravity black 
hole solutions |l], ||, ^ have been understood microscopically to be described by D-branes 
wrapping supersymmetric cycles 0, M. Nonetheless, all microscopic BPS configurations 
cannot correspond to single centre solutions in the effective N = 2 supergravity description, 
and it has been understood @ that the latter then describe composite bound states of BPS 
black holes @) Q §] • The detailed description of these solutions and their domain of stability 
in moduli space was very important in order to correctly reproduce the corresponding 



microscopic results P, 10 1. 



In the aim of generalising this understanding to more realistic non-supersymmetric 
black holes, the first non-trivial step is to consider non-supersymmetric extremal black 
holes. One can classify extremal black hole solutions into two main classes: the over- 
rotating class, for which the angular momentum saturates the extremality bound, and 
the under-rotating class, for which the electro-magnetic charges saturate the extremality 
bound, and which includes the supersymmetric solutions. Given that a stationary space- 
time defines a time fibration over a space-like three-dimensional base, one may distinguish 
the various classes by their corresponding base space. The characteristic of all known under- 
rotating extremal solutions is that the three-dimensional base space is flat, i.e. M^ with 
the Euclidean metric, whereas the over-rotating extremal solutions all admit the singular 
three-dimensional base of the extremal Kerr solution. In this paper we will only discuss 
under-rotating solutions admitting a flat three-dimensional base space. ^ 

Within the class of under-rotating solutions with an M^ base, there are two known 



interacting systems of non-BPS black holes. The composite non-BPS system [11, 12, 13 1 



describes the interactions of black holes that are non-BPS in isolation, while the almost- 



BPS system |14, 15, 16], allows for configurations of centres that are be both BPS and 
non-BPS in isolation. Both these systems have been studied in some detail and, in several 
cases, the known explicit solutions are general enough to allow for the most general solution 
to be obtained by dualities [i.e. symmetries of the three-dimensional theory). 

An important feature of the existing formulations of both these non-BPS systems is 
the existence of distinguished directions in charge space, such that some are associated to 
harmonic functions as in the BPS system, whereas others are associated to functions solving 
Poisson equations. The different directions being then intrinsically inequivalent, one cannot 
straightforwardly rotate one into another by duality, as for BPS solutions. It is therefore 
customary to solve these equations and construct solutions parametrised by integration 
constants that can only be related to charges and asymptotic moduli a posteriori. Of 
course, one may straightforwardly apply a duality rotation on the solutions to change the 
identification of pre-determined charges, thus covering all possible charge configurations. 



^In the five dimensional uplift these solutions feature more elaborate base spaces, which however still 
include a flat M? . 



This process is however not only cumbersome, but more importantly it obscures the overall 
structure of these systems, which should be manifest in a fully covariant formulation. 

A further motivation for generalising the existing formulations is that non-supersym- 
metric under-rotating extremal black holes seem to admit a microscopic description similar 
to their BPS cousins [17, |l^]. Understanding the domains of existence and stability of 



these solutions will eventually be important in order to understand in more detail their 
microscopic description. However, despite the existence of explicit solutions describing 
bound states of such non-BPS black holes, their domain of existence in moduli space has 
not been studied in detail. One of the main obstacles in carrying out this program originates 
from the property that these solutions are described in terms of parameters that are not 
the physical charges themselves, as explained above. In order to solve this problem, it is 
important to obtain the general solutions associated to fixed charge configurations. Being 
able to do this in a straightforward manner requires the definition of covariant equations 
that are not constrained to a fixed duality frame. 

In this paper, we provide a manifestly duality covariant formulation of both the com- 
posite non-BPS and the almost-BPS systems that unifies their description, by rewriting 
these systems of linear differential equations in terms of duality covariant quantities, similar 
to the BPS system ||6|. Although these systems of equations turn out to be more compli- 
cated than the BPS one, the formalism permits to compute the most general solutions 
associated to fixed charge configurations. The construction of explicit solutions within 
these systems will be discussed in a forthcoming publication. 

Somewhat surprisingly, we find that both systems are naturally formulated by consid- 
ering space dependent translations along abelian isometries of the scalar manifold. As all 
equations are written in terms of objects transforming linearly under electromagnetic dual- 
ities, it is natural to introduce local generators for the abelian subgroups corresponding to 
these isometries. The relevant equations in both the composite non-BPS and the almost- 
BPS can then be expressed in terms of a covariant derivative that contains a nontrivial 
connection in the Lie algebra of the relevant abelian isometries. 

In a specific frame, these abelian isometries coincide with the isometries that generalise 
the action of T-dualities combined with large gauge transformations of the p-form gauge 
fields present in string theory, for continuous parameters. In general, they will define 
abelian subgroups that are homomorphic to the latter by similarity transformations. For 
simplicity, we will refer to them as T-dualities in this paper, despite the fact that they are 
not in general associated to any duality in string theory. The set of all abelian isometries 
in the scalar target space of symmetric models has been considered before in the context 



of (non-)BPS black holes in different guises, see e.g. [^, 20]. 

The paper is organised as follows. The remainder of this introductory section is devoted 
to an informal presentation of the method we use to obtain the duality covariant formulation 
of the composite non-BPS and almost-BPS systems. In doing so, we also explain in more 
detail way how the T-dualities arise in both systems and in fact allow us to define their 
action on the relevant quantities without introducing explicit matrix representations. In 
section |2| we present the general action of arbitrary T-dualities in terms of symplectic 
vectors and discuss the corresponding decomposition of the charge vector space. In section 



m we consider the three dimensional Euchdean coset non-hnear sigma model describing 
stationary black hole solutions. We then explicitly solve the nilpotency conditions on the 
scalar momentum as a Lie algebra element for both the composite non-BPS and the almost- 
BPS system, to obtain the first order flow equations describing each class. We then go on in 
sections ^ and ^ to rewrite these flow equations as a linear system of Poisson equations for 
a set of symplectic vectors parametrising the solutions in an arbitrary duality frame. The 



reader interested in applications can find summaries of the two systems in sections 4.5 and 



5.5 , which are self contained and only require the definition of T-dualities in section 2.2 
We conclude in section y, where we discuss some of the implications of the local T-dualities 
acting on both systems and comment on further generalisations. 

1.1 Overvievif of results 

In order to find systems describing stationary black hole solutions, we consider the reduction 
along time to a three-dimensional Euclidean theory. In this setting, one can straightfor- 
wardly dualise all vector fields to scalars, to obtain a non-linear sigma model coupled to 
Euclidean gravity. In this paper, we only consider J\f = 2 supergravity theories with a sym- 
metric scalar space, so that the scalar fields of the resulting non-linear sigma model over a 
pseudo-Riemannian symmetric space. In the case of solutions with a flat base space, there 
is a powerful method for obtaining first order fiow equations that solve the full equations 
of motion, starting from the observation that the Einstein equation 

Tr(P^P,) =i?^, = 0, (1.1) 

can be solved by assuming the scalar momentum, P, to be a nilpotent element of the Lie 
algebra. It appears then that the scalar fields equations of motion reduce to a solvable 
system of differential equations. Based on standard group theoretical considerations, one 
shows that the nilpotency of P implies that there exists a second element of the algebra, 
h, such that an eigenvalue equation of the type 



n(h-z)P = 0, (1.2) 



i=l 



holds, where hP = [h, P] and n is a positive integer. It then follows that P must be a linear 
combination of eigenvectors of h with positive eigenvalues. In section ^ we consider the 
relevant eigenvalue equations for both the composite non-BPS and almost-BPS systems, 
giving the precise solutions for both the momentum P and the auxiliary fields h. For the 



composite non-BPS system, this analysis has been given in [13|, but is included here for 
completeness. The corresponding flow equations for the almost-BPS system have been 
derived for the STU model in [^ by reduction of the second order equations of motion, 
while we present the general analysis for any symmetric model. 

The presence of the aforementioned auxiliary fields is a central feature of both systems, 
as will become clear from our treatment. In fact, we find it most convenient to keep them 
throughout, as variables characterising the solutions, despite the fact that one may in 
principle solve for them in terms of physical fields. In the systems we shall consider, h can 



be parametrised in terms of a distinguished direction in the charge vector space (plus a 
phase for the almost-BPS system). This distinguished direction is associated to a constant 
so-called very small vector, or in other words a one-charge vector, in the sense that it can 
always be brought by dualities to a canonical form where it has only one charge, as a pure 
D6 charge for example. However, one must keep this very small vector arbitrary in order 
to be able to consider generic charge configurations. 

Treating the system in terms of electromagnetic charges and auxiliary vectors implies 
that one has to work with real quantities, rather than the complex quantities that appear 
naturally in the Euclidean 3-dimensional non-linear sigma model. It turns out that this 
is essential in solving the system in terms of local functions, since a linear structure only 
appears after writing all flow equations in the real basis. As mentioned above, this change 
of basis, that permits to solve non- linear first order equations in terms a solvable system 
of linear differential equations, is the main technical result of this paper. 

In order to appreciate the importance of the real formulation, the simple example of 
the analogous situation for the BPS system is instructive. In this case the auxiliary element 
h in ( |1.2| ) is much simpler, as it is parametrised by a single phase, e*°. Taking the static 
case for simplicity, (1.2) leads to the standard BPS equations 

d,,U = r-2e^Re(e-^"Z(r)) , drf = r-"^ e^ e'"" Z\T) . (1.3) 

Here, r is the distance from the horizon, e is the single function parametrising the metric, 
t* are the vector multiplet scalars and Z(T), Zi(^) &i^6 the central charge of the charge 
vector r, and its Kahler derivative. In this form, it seems a rather non-trivial task to 
solve (|1.3| ) explicitly. It is however possible to combine the scalar degrees of freedom in a 
single vector, e~^e~*"V, where V is the so called symplectic section, which parametrises 
the scalars i*. The two equations (|1.3D can then be written as in [o]. 



2 9Jm (e-^e-^°V) = ^ , (1.4) 

using standard special geometry identities. It is now trivial to solve the last equation as 

2Im(e-^e-*"V) = --H, (1.5) 

where T-L \s a vector of harmonic functions, whose poles are identified with the charges T. 
One still has to solve a set of algebraic equations to obtain the scalars i* and e^ in terms 



of harmonic functions H, but it is important to stress that (1.5) is equivalent to a solution 
for the physical scalars. 

The analogous computation for the non-BPS case includes the auxiliary vector para- 
metrising the element h and is therefore considerably more involved. For the restricted 
case of single centre under-rotating black holes, this was done in detail in [^] and the 



result was given in terms of the objects appearing in (L5), together with the vectors in 
( |1.2| ) above. This class is naturally part of both the multi-centre systems considered in 
this paper and its description is crucial for understanding the structure of the more general 



classes. Explicitly, for the single centre class we have |22, 13| 



2Im(e-^e— V)=-?^ + 2ff|i? + ^i2% (1.6) 



where Ti is again a vector of harmonic functions describing the charges of the black hole, 
while three new objects appear, namely the function M and the two constant very small 
symplectic (pseudo-charge) vectors R and R*, which are not mutually local, i.e. {R,R*) 7^ 
0. 

The presence of two different very small vectors in ( |1.6D can be somehow surprising, 
but it is important to note that they do not define independent parameters. Indeed, we 
find that these vectors are determined by the electromagnetic charges T, up to duality 
transformations leaving T invariant [0| , so that they are determined in terms of T and the 
asymptotic scalar fields t^. Conversely, one can view the charges, F, and the harmonic 
vector Ti as being constrained to lie in a Lagrangian subspace determined by R and R*. 
The terms proportional to the two very small vectors are worth discussing in more detail. 
First, note that the term proportional to R is simply a projection on that component of 
Ti with an additional factor of 2, implying that this particular component appears with a 
flipped sign. This is a very general feature of non-BPS solutions that has been observed 



in many examples in the literature [^, 24, |2^. On the other hand, the component along 
R* in (pB) contains the function, M, which represents the only genuinely new term in the 
expression for the scalars, and is constrained to be a dipole harmonic function characterising 



the angular momentum of the under-rotating non-BPS black hole |15]. 

In sections |^ and |5|, we explicitly solve the composite non-BPS and almost-BPS non- 
linear first order systems derived in section y in terms of real vectors of local functions. As 
we shall see, the generalisation to multi-centre systems can be performed by allowing one 
of R or R* to vary in space, while keeping their symplectic product fixed {R, R*) = 4. The 
space dependence of the non-constant vector can then be reabsorbed into a space dependent 
duality transformation that leaves the constant very small vector invariant. Because this 
duality transformation lies in an abelian subgroup that is conjugate to the group of T- 
dualities by similarity transformations, we shall simply refer to them as T-dualities. In 
particular, the composite non-BPS system is described by the T-dualities, T+, defined 
as leaving R invariant, while the almost-BPS system is described by the T-dualities, T~, 
leaving R* invariant. One then shows that the symplectic section V still takes the form 
( |1.6| ), with Ti constrained to lie in the same Lagrangian subspace determined by R and 
R* . Since one of the very small vectors is not constant, Ti is not harmonic anymore, but 
satisfies the Poisson equation 

d • d{exp[-T'^]n) = dT^ A *dT± exp[-T±]?^ , d * dT± = , (1.7) 

where the second equation imposes that the n^ parameters of the T-dualities are given 
by arbitrary harmonic functions. The function M in (|1.6D is also specified by a T-duality 
covariant equation, given by 

i^diu-dM={'H,dn-2dT^n), (1.8) 

which also fixes the angular momentum one- form, lo. Although it is not manifest from 
these general equations, the graded structure of the vector space is such that the system is 
solvable. Because the two T-dualities in the right hand side of ( |1.7D have a nontrivial kernel. 



it follows that the source in ( |1.7D does not contain some of the components of exp[— T ]?^, 
which therefore includes both harmonic and non-harmonic components. It turns out that 
T^ act as raising or lowering operators, so that only the harmonic components turn out 
to source the non-harmonic ones, i.e. the non-harmonic components of exp[— T^]?^ do not 
source themselves. These properties will be discussed in detail in sections ^ and |5|. 

By definition, T-dualities are represented in terms of symplectic matrices in the 2(n^ + 
l)-dimensional vector space of charges. However, as will be shown in detail in the following 
sections, the n^ parameters of a T-duality can be arranged into a symplectic vector obeying 
a number of constraints. One can therefore write the action of the corresponding T-duality 
in terms of this vector, together with R and R*, using the symplectic product and the 
quartic symmetric invariant I4 (which defines the entropy of extremal static black holes). 
In this way, we obtain explicit expressions for the sources in (1.7), in terms of exp[— T ]?^ 
and the parameters of the T-dualities, that can be evaluated explicitly once a model and its 
associated quartic invariant are specified. These equations reduce to the systems introduced 
in ||ll|, 0], for a particular choice for the two very small vectors above. 

This concludes our short presentation of the main results in this paper. For the con- 
venience of the reader, we provide an account of the results appearing in the following 
sections, which can be read independently of each other with the exception of section |2|, 
that is basic to most applications. In section |2| we give a detailed discussion of the general 
T-dualities. We show that given two very small vectors R and R* that do not mutually 
commute, one can explicitly define a graded decomposition of the symplectic vector space, 
on which the T-dualities T^ can be respectively defined as raising and lowering operators. 
The structure of this decomposition is essential for all applications in this paper. 

We then go on in section ^ to define the two systems of non-BPS multi-centre black 
hole solutions in terms of the non-linear sigma model in three dimensions, obtained after 
dimensional reduction over the time direction. After an overview of the main properties 
of this three-dimensional Euclidean theory, we discuss in detail the eigenvalue equations 
( |1.2| ) for the two systems at hand. This results to two sets of first order flow equations that 
completely describe the composite non-BPS and almost-BPS systems in four dimensions. 

Sections ^ and || are mirror copies of each other, wherein we present in detail the change 
of variables that transforms the flow equations of section ^ to two linear systems. The 
reader can find a concise summary of the two multi-centre systems in the real formulation 



in sections 4.5 and p^ respectively, where we also discuss some general properties of the 



solutions, deferring a more detailed presentation and explicit examples for a forthcoming 



publication ||26(| . 



2. T-dualities 

In this section, we provide a detailed discussion of the abelian isometrics on the scalar 
target space of A/" = 2 supergravity,^ that leave a given charge vector invariant. As will 
be shown in later sections of this paper, the precise action of these isometrics, in their 



■^We refer to the appendix for a concise review of our conventions and notations. 



most general form, is of central importance in the construction of multi-centre non-BPS 
black holes. We start with an informal discussion of the simplest example of such abelian 
isometrics. Subsequently, we review some properties of symmetric special Kahler spaces 



with a cubic prepotential in section 2.1. Section 2.2 is devoted to the definition of the 



abelian isometries we shall refer to as T-dualities and their explicit action in terms of the 



associated real vectors. In section 2.3 we discuss the realisation of the same isometries in 
the complex basis defined by the central charge and its Kahler derivatives. 

Throughout this paper, we study extremal multi-centre black hole solutions in A/" = 2 
supergravity coupled to n^ vector multiplets labeled by an index i = l,...,n^, whose 
scalar fields, t*, parametrise a symmetric special Kahler target space, A^4. These spaces 
were classified some time ago |2^ and include minimally coupled vector multiplets, which 
are not of interest in this work, and theories with a cubic prepotential, specified by a 
completely symmetric tensor cijk [^] (cf. ( |A.5D ). The target space geometry is governed 
by a Kahler potential (cf. ( |A.7| )), which manifestly depends only on the imaginary part 
of the scalars, t*. It follows that the real parts of the scalars are coordinates along n^ 
isometries of the scalar manifold, acting as 

f^f + k\ (2.1) 

where A;* is a vector of n„ constant real parameters. These isometries are generic in all cubic 
models and clearly form an abelian algebra. For theories originating from Calabi-Yau string 
compactifications, the operation (2.1) can be viewed as large gauge transformations on the 



higher dimensional tensor gauge fields along internal cycles, combined with T-dualities. 

While this description is useful in characterising the symmetries themselves, for reasons 
that will become clear below, in this work we are interested in the embedding of these 
isometries in the symplectic group, which acts on the electric and magnetic gauge fields 
in four dimensions. The convenient variable to use in order to make the action of the 
isometries in ( |2.lD transparent is the so called scalar symplectic section, V, which is a 
somewhat redundant way of repackaging the scalars, as 

v=r; 1 , t' = |^, (2.2) 

where Fi are the derivatives of the prepotential with respect to the X . The index 
/ = {0, i}, runs over one more entry than n« and enumerates all the gauge fields in the 
theory, i.e. the vector multiplet gauge fields and the graviphoton. Note that V changes 
under Kahler transformations by a phase and is subject to the constraint (A. 6), so that 




it encompasses only 2n„ degrees of freedom, identified with the physical scalars t*s. The 
advantage of this variable is that, unlike the physical scalars, it transforms linearly under 
electric/magnetic duality transformations, in exactly the same way as the electromagnetic 
charges. 

For instance, the isometries in ( p.l| ) are described by a linear transformation acting on 



the charges as 



exp[T;; 



/pO\ 
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p' 






-^ 


Qi 




\qo) 


\ 



Qi + Cijkk^p^ + I Cijkk^k^p^ 
\qo - k'qi - i Cijkk^k'^p' - Icijkk'k^k'^p^ J 



(2.3) 



where we defined the abehan generators T^ for later convenience. One can now easily 
verify that the same operation ( p.3| ) acting on the section in ( |2.2| ) leads to ( p.lj ) for the 
physical scalars. In this formulation, the connection of the isometries ( |2.1| ) to higher 
dimensional gauge transformations is more transparent, since it has a natural action on 
the electromagnetic gauge fields. Moreover, this particular set of abelian transformations 
also arises in the form of spectral flows in conformal field theories describing black holes 
microscopically. 

The crucial feature of the symplectic embedding of the abelian isometries is, however, 
that one may generate an infinite number of inequivalent sets of abelian isometries by 



conjugating the matrix exp[Tfc] in (2.3) by a general [/-duality transformation, as in pO| . 
^ These sets of isometries are more complicated than the one in ( |2.1| ) and do not commute 
with it. From a higher dimensional point of view, some of these more general isometries 
can also be viewed either as large gauge transformations conjugated with generic T- and/or 
[/-dualities, or as (generalised) spectral flows in a dual conformal field theory. Here, we 
refer to them simply as T-dualities for brevity and we focus on the case of symmetric scalar 
manifolds, which allows for the most general transformations to be described explicitly. 

The representation of a generic T-duality in terms of matrices is of course a rather 
tedious task, which can be circumvented in a natural way, intrinsically tied to the systems 
of non-BPS black holes we consider. The crucial observation is that there is always a graded 
decomposition of the vector space in four components, generalising the clear distinction 



between the various components in (2.3), based on their transformation rule under T- 
dualities. Indeed, general T-dualities act consistently on each component of the charge 
space with a fixed homogeneity in the parameters /c*, which can never exceed three. In 
particular, there is a distinguished direction that is invariant under the action of any given 
T-duality, as for example the electric charge go is left invariant in ( [2.31 ). The qo charge is 
rather special, since it is an example of a so-called very small vector. We will recall the 
precise definition of such a vector in what follows, but loosely speaking a very small vector 
can be defined as a 'one charge vector', in the sense that it is [/-dual to a pure qo charge. 
Clearly, the distinguished direction that is left invariant under a generic T-duality must 
then always be a very small vector, given that all such transformations are [/-dual to the 
above example. 

The relevance of very small vectors for extremal non-BPS solutions arises already in 
the single centre class [22, 13, ^], which is naturally described in terms of an auxiliary 
pair of mutually nonlocal constant very small vectors, constrained by the physical charge 



■^Following a common misuse of language, we call [/-duality transformations all the continuous isometries 
of the scalar field symmetric space, independently of whether they define or not actual string theory dualities 
for appropriate integral coefficients. 



vector. As we will show explicitly in sections ^ and g below, the corresponding multi- 
centre systems are naturally described by promoting one of these very small vectors to be 
not constant. In order to see the connection to T-dualities, consider the very small vector 
R defined such that its only non-vanishing component is qo = 4. A general very small 
vector can be parametrised as 

S = c(l,s'; ^Cijksh'', -Icijks's^s''^ , (2.4) 

where c and s* are allowed to take singular values as long as the components of S are well 
defined in the limit. '^ It follows that a general very small vector R* satisfying {R,R*) = 4 
can be parametrised as 

R* = exp[T+] R* = (l, s' ; ^c.jkS^s'^, -^Q.-fc/s^s'^)^ , (2.5) 

where the only non- vanishing component of Rq is p^ = 1. In the case of a constant vector 
R*, (|2.5| ) is simply a convenient parametrisation, but in the more general case when R* is 
not constant, one can assume the parameters s* to be functions of space, to obtain 

dR* = dT+ R* = T+^ R* , (2.6) 

where we used the abelian property of ( |2.3D . As we will show explicitly in later sections, 
the composite non-BPS system is naturally characterised by two very small vectors, one 
constant R and one non-constant R* , which have a non- vanishing constant symplectic 
product {R,R*) = 4. The almost-BPS system is similarly described by two very small 
vectors, only the role of R and R* are interchanged. The constant very small vector is left 
invariant by the relevant T-dualities (as for example (^]^)), which are therefore different 
for each system. The non-constant very small vector can be expressed as in ( |2.5D for a 
constant vector Rq in the composite non-BPS system (or respectively Rq for the almost- 
BPS system), and it follows that it satisfies (|2.6| ). In the rest of the paper R and R* will 
be generic very small vectors, and the associated T-dualities will define abelian subgroups 
conjugate to the one described in 



2.1 Symmetric special Kahler spaces 

In this paper we consider Af = 2 supergravity theories defined in [^] for which the special 
Kahler target space, A^4, is a symmetric space and that can be obtained as Kaluza-Klein 
reductions of corresponding five dimensional theories.^ In this case, Ai^ is a coset space of 
the four-dimensional duality group, G4, by its maximal compact subgroup C/(l) x K4 

M4 = (C/(l) X i^4)\G4 . (2.7) 

For the class of theories we consider, the scalar target space is a symmetric space even 
after dimensional reduction/oxidation to three/five dimensions, so that ( p. 71 ) is part of the 



*For example one recovers R in the limit c ^- 0, s' — )■ 00 with —c^CijkS^s-' s'' — 4. 

^This excludes theories with minimally coupled vector multiplets, which do not include systems of the 
type we consider here. 
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sequence of embeddings 

K5\G5 -^ (f/(l) X K4)\G4 -^ {SL{2) X G4)\G3 , (2.8) 

where by Gd, Kd we denote the duahty group and (part of) the isotropy group in d 
dimensions respectively. Note that the divisor group in three dimensions is non-compact 
because we consider the time-hke reduction to three dimensions,^ as that is relevant for 
the applications we consider later on. Note that K4, is the compact real form of G5, by 
property of A/" = 2 super symmetry. 

One can always define a set of vielbeine associated to the Kahler metric gij on M.4^ 

f^ = ele'^\ (2.9) 

such that the constant symmetric tensor 

Cafec = ie'^4efee^Cijfc, (2.10) 



where Cjjfc is the G5 invariant tensor defining the prepotential (cf. ( |A.5 )), is left invariant by 
Ka,. Then, the contravariant symmetric tensor c°^'^ in the conjugate representation satisfies 



the Jordan identity [28 



Cf(abCcd)gC''^^ = ^^f^Cbcd) ■ (2-11) 

In a complex basis, the Lie algebra of G4, denoted 94 (and respectively €4 for K4), 
naturally decomposes as 

04 = u(i) e «4 e c"" . (2.12) 

It follows that the relevant parameters are given by those corresponding to the elements of 
^4, denoted by G°fe, a real scalar 7 and a complex vector A^. The corresponding algebra is 
realised in terms of anticommuting parameters with the nilpotent differential^ 

5Aa = -G\Ab + 2i7A, ,57 = ^A'^A^ , 

6G\ = G\G% + c'"''cMe^,A''+A'^A, + ^A^A,<5,« . (2.13) 

Note that the statement of invariance of the tensor Cabc under i^4 implies that ^[2] = 
Icabc^"' Z^ Z'^ is K4 invariant for any vector 2°" transforming in the relevant n„-dimensional 
complex representation of i^4. One can check that the variation of G°f, in ( 2.13| ) indeed 



leaves invariant the cubic norm J\f[Z] = M\Z\ for an anticommuting A^- 

The invariance of the cubic norm M\Z\ can be used to define duality invariants and 
restricted charge vectors, a concept that is of central importance for the applications we 
consider later in this paper. First, we introduce the quartic invariant for a charge vector 
r, in terms of its central charges, Z = Z{T), Za = Za{T), as 

/4(r) = (ZZ -ZaZ^^f -CeabZ''Z''c^^'^ZcZd + 4ZM[Z]+4ZM[Z]. (2.14) 



"For a spacelike reduction, one obtains SU(2) x K3, where Ki is the compact real form of G4. 
^The nilpotent differential 5 in these equations acts as SXaT — ^ XaXb[T ,T^] on a Lie algebra 
element expressed in a basis of generators T using a set of anticommuting parameters Xa ■ 
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This expression can be verified to be invariant under the 04 generators of ( 2.13| ), so that it 



is moduli independent. This is manifest by the corresponding real form of this invariant, 
which is given solely in terms of charges by 

2 2 

= -{p°qo + P\if + 3 go CijkP'p>P^ - - P° c'^^qimk + CijkP'p'' c'^^'qiqm , (2.15) 

where we also defined the completely symmetric tensor t Q for later reference. Again, 
one can easily check that ( p. 15 ) is invariant under the example T-duality in (^]^), and it 



is more generally invariant under an arbitrary G4 transformation. 

We are now in a position to introduce the concept of charge vectors of restricted rank. 
A generic vector leads to a nonvanishing invariant ( 2.14 )-( p.l5D and is also referred to as a 



rank-four vector, due to the quartic nature of the invariant. Similarly, a rank-three vector, 
Fj, is a vector for which the quartic invariant vanishes, but not its derivative. An obvious 
example is a vector with only p* 7^ and all other charges vanishing, so that the derivative 
14 (Fj) is nonzero and proportional to the cubic term Af\p]. 

There are two more classes of restricted vectors, defined analogously as rank-two 
(small) and rank-one (very small) vectors. A rank- two vector, T2, is defined such that 
both /4(F2) = 14 (F2) = 0, and a simple example is provided by a vector with all entries 
vanishing except the p*, with the additional constraint that Af[p] = 0. Finally, a very small 
vector, Fj, is defined such that 

/4(F,) = I^(F,) = 0, 
i/4(F„F„F,F) ^ it*^^^QF,MF,^FpFQ = -(F„F)2, (2.16) 

for any vector F. In this case, we can also give a general definition in terms of the complex 
basis, which is in fact independent of the values of the scalar fields. In this paper we will 
often make use of a rank one vector, R, that we choose without loss of generality such that 
|Z(i?)| = 1. One shows that such a very small vector satisfies 

ZiR)=M[n], ZaiR) = na, (2.17) 

where J\f[0,] is a phase by construction. The remaining central charges Oq are such that 

-c°''"'nbnc = J^[n]na, f)"f^a = 3. (2.18) 

A general very small vector can be obtained by rescaling both 7V[f]] and fla by a real 
function. Examples of very small vectors were already given above, as vectors where only 
the qo or p^ component is nonzero, while the parametrisation given in ( p. 5]) is generic up 
to a possibly singular rescaling. 

2.2 Freudenthal ternary algebra realisation of G4 

We now proceed in describing the duality group G4, as defined above, in terms of real 
vector parameters. This is essential for discussing the T-dualities, which are contained in 



12 



G4 as subgroups and can therefore also be described in terms of real vector parameters in 
the general case, similar to the example ( |2.3[ ) above. 

The central object for the definition of G4 in the real basis is the quartic invariant 
in ( 2.15| ) and its derivatives. It is convenient to define a symplectic vector out the first 
derivative, /4(r), of the quartic invariant so that 



(r, i^(r)) = 4/4(r) , /l(r, r, r) = 6/^(r) 



(2.19) 



Using the definition ( 2.15 ) of the quartic invariant and the properties of the rank-three 
symmetric tensor Cijk, one shows the following quintic identity 



/i(r,r,i^(r)) = -8i4(r)r, 



(2.20) 



for a generic charge vector T. This identity is equivalent to the property that the Freuden- 
thal ternary product 



iX,Y,Z) ^hi{X,Y,Z) + ^{X,Y)Z -^{Z,X)Y + ^{Y,Z)X 



(2.21) 



satisfies the four axioms defined in pO], and inversely, any Preudenthal ternary product is 
necessarily of the form (2.21), for a completely symmetric rank four tensor satisfying ( 2.2C| ). 
One can therefore define the 34 Lie algebra as in ISH]. We shall not use the Freudenthal 
ternary product, but rather the quintic identity ( 2.2Cl| ). We refer to [^] for the more formal 
definition of the 94 Lie algebra from the ternary product itself. It is straightforward to 



combine ( 2.20 ) with the symmetry properties of the sextic invariant {I'{Ti),I'{T2)) to show 
that for any two vectors Ji and J2, the linear transformation 



g(Ji, J2)r = i/^Ji, J2,r) - Ji(J2,r) - J2(Ji,r) , 



(2.22) 



preserves both the symplectic product and the quartic invariant. One concludes that ( 2.22| ) 
defines a generator of 04, and all 94 generators can in fact be defined in this way. It follows 
that the Lie algebra takes the form 



[giJl,J2), g(J3, J4)] = g(g(Jl, J2) J3, J4) + g(J3,g(Jl, ^2)^4) 



(2.23) 



as shown in [pQ] . 

A special case arises for two rank 1 vectors, denoted R and R*, which are assumed to 
be mutually non-commuting. In this case, one can define the corresponding 34 generator 
as in ( ^321) 

hrV = {R,R*)-^(h'^{R,R*,T) + {T,R*)R- R*{R,r)) , (2.24) 



which is central in the description of T-dualities. It is clear from ( p. 16 ) that for any rank 
1 vector R (or respectively R*) one has 



I'4{R,R,T) =4{R,T)R . 



(2.25) 
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This generator admits therefore R and R* as eigenvectors, with eigenvalues +3 and —3 
respectively, and the remaining eigenvectors of Iht can be characterised as follows. Using 
( |2.20 ) and ( 2.16| ), one can show that 



jli{R,R*,li{R,R*,T)) = {R,R*fT + 3{R,R*){{T,R*)R + R*{R,T)) , (2.26) 



(2.26) from which follows the action of the square of hj-, as 

h^r = r + 8{R, R*)-\{r, r*)r + r*{r, r)) . (2.27) 

This equation implies that R and R* are the unique eigenvectors with eigenvalues +3 and 
—3 respectively, and also leads to the characteristic equation 

(h^ - lOh^ + 9)r = (hT - 3)(hT - l)(hr + l)(hT + 3)r = . (2.28) 

In view of the fact that h-p is symplectic, it follows from ( |2.27] ) , ( p.28D that the 2n„ + 2 
electromagnetic charge vector space decomposes into 

j^2n„+2 ^ j^(-3) ^ (Jj^n«)(-i) q (^^^ )W ^ j^(3) ^ (2.29) 

where the two distinguished vectors R and R* are by definition the components of grade 3 
and —3, respectively. This decomposition is clearly relevant to the T-dualities as described 



in ( p. 290 , as it allows to identify four eigenspaces, based on two very small vectors. The 
Hy eigenvectors of eigenvalue +1 and the n^, eigenvectors of eigenvalue —1 can be obtained 
by defining the corresponding projectors to the four eigenspaces of hy, as hj-r*"' = nF'"^ 
for n = —3, —1, 1, 3, i.e. 

r(^) = {R,R*)-\T,R*)R, 

rW = ^r + ^{R,R*)-^(^IiiR,R*,T)-3{T,R*)R + R*{R,r) 

r(-^' = lr-l-{R,R*)-^(ll'^{R,R*,T)-{T,R*)R + 3R*{R,r)' 

T^-^^ = {R,R*)-^{R,T)R* . (2.30) 

Note that the r*='=^' can simply be identified as the solutions to 

h'^{R, ii*, r(±^') = ±{R, i?*)r±i> . (2.31) 

These expressions will be very useful in evaluating the action of T-dualities on general 
symplectic vectors in the following sections. A further practical advantage of this decom- 
position is the fact that all inner products must respect the grading, leading to strong 
constraints on the possible nontrivial combinations. For instance, the grading implies that 

I'{r^-'\r'--'\R*) = , /'(FW,F(i',i?) = , (2.32) 

since there is no vector of weight ±5 that these cubic terms could be equal to. Similar 
considerations apply to scalar products, which necessarily vanish unless the sum of grades 
of the vectors involved vanishes. 
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In addition to the decomposition ( 2.29| ) of the vector space, the generator h-p imphes a 



corresponding decomposition of the duahty group generators. Indeed, hy commutes with 
05 C 04 and defines the fohowing graded decomposition of 04 

04 ^ (R"")(-^) e (0[i 05)'°' e (]R"")''> , (2.33) 

where the qIi corresponds to h^ itself. Clearly, the 2?!^, generators of eigenvalue ±2 with 
respect to hy can be used as raising and lowering operators on the eigenspaces in ( 2.2g| ). 



As the reader might already understand, these grade 2 generators are related to the trans- 
formations ( p.3[ ) by similarity transformations in G4. 

In terms of the explicit expression ( 2.22| ) for the action of 04, one may consider any 



grade —1 vector of parameters A;'"^' to define the grade 2 generators as 

T+r ^-^iR, R*)-' (lliiR, k^-'\T) + (r, A;'-^')^ - k^-'' {R, r)) • (2.34) 

This generator is manifestly of grade 2 because of the grading of R and fc*'^' themselves 



and the algebra ( 2.23 ). It is convenient to write it in a way that makes the grading explicit 



T+r = {R,R*)-^(k^-''>{R,T'^-''>) - ^li{R,k'--''>,T^-''>) - (r<^\ /t'-^O^) , (2-35) 



where we used the projections in ( |2.3[1| ) and the fact that /c*"^' is of grade (—1). All these 
generators clearly commute between themselves for different A;'"^''s. Similarly, one defines 
the grade —2 generator in terms of a grade 1 vector A*^' 

T^r = ^{R,R*)-'(^^liiR\k^'\T)-{k^'\r)R* + k^'\T,R*) 

?*\-i/^i-(i)/r(3) p*\_uir'^p* i.(i) rW\-/i.w r(-i)\ 



= {R, R*)-^ (k'-'^T'^''\ R*) + -4(i?*,A«,r«) - (A;(^',r(-^')i?*j . (2.36) 
The normalisations we have chosen are such that 

T+R* = k'--^'' , T^R = k^'^ , (2.37) 

while one easily computes that 

T+R = , T^R* = . (2.38) 

In this form, one easily computes that these generators are nilpotent of order 4, as 

(T±)^r = , (2.39) 



consistent with the grading (2.29), which only allows for four eigenspaces. Explicitly, we 



find the following expressions for the two sets of generators 

(T+)2r = -^{R,R*)-^(^l',{R,k^-'\k^-'^){R,r)+h{R,k^-'\k^-'\r)R^ , (2.40) 
(T+fV = -^{R,R*)-''hiR,k^-'\k^-'\k^-'^){R,T)R , (2.41) 



(T^)2r = -{R,R*y'^(^l'^{R*,k^'\k^'^){r,R*)-h{R*,k^'\k''\T)R*j , (2.42) 
(T^)3r = -^{R,R*)-^h{R*,k^'\k^'\k'^'^){T,R*)R* , (2.43) 
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and moreover 

[hr,T±] = ±2T±, (2.44) 

as in ( p.33 ). 



Finally, one also computes using ( p. 23 ) and the grading that for any grade —1 vector 
e and grade 1 vector / 

[T^ TT] = 1—1— f g(e, /) + (e, mr) . (2.45) 



^f' 2{R,R*) . 
One can straightforwardly check that 

'g{e,f) + UeJ)hT)R = 0, (g{eJ) + UeJ)hT)R* = 0, (2.46) 



so that the latter transformation lies in the 05 C 04 subalgebra, consistently with the 
graded decomposition ( ^.33]) . One can indeed check that these transformations preserve 
the cubic norm /4(i?, r'"^',r*'^\r'"^') for an arbitrary grade —1 vector r'"^^ It turns out 



that the identity ( 2.20 ) implies the associated Jordan identity 



64 

liiR*,liiR,p-'\r^-'^),liiR,T^-'\T^-'^)) = —{R,R*)h{R,T^-'\T^-'\T^-'^)T^-'^ , (2.47) 

O 

which generalises (|2.11| ). These equations are clearly valid upon replacing R with R* and 
r<"^' by a grade (+1) vector F'^* throughout. 

One may now use the above formulae to identify T^ with T-dualities explicitly. In- 
deed, one can easily check that upon identifying R with the very small vector whose only 
nonvanishing component go and R* with its magnetic dual along p^, the exponentiated 
transformations 

exp[T+] = 1 + T+ + 1 (T+)2 + 1 (T+)3 , (2.48) 

are identical to the spectral flow shown in (|2.3| ). The corresponding set of generators T" 
then generate the T-dualities one obtains by conjugating ( |2.3| ) by an electric/magnetic 
duality and leave R* invariant. The grade (—1) and (+1) components are then easily seen 
to be given by the magnetic, p*, and electric components, Qi, respectively. 

In the general case, we can identify all possible sets of T-dualities as given by a choice 
of Rot R*, as above, as the generators T^ are entirely determined by the rank 1 vector they 
leave invariant. Indeed, the characteristic feature of these abelian subgroups is that there 
is always a unique (up to rescaling) very small vector {e.g. R) that they leave invariant, 
whereas they act transitively (up to a rescaling) on the set of very small vectors {e.g. R*) 
that are not mutually commuting with the former. In the specific example of ( |2.3| ), any 
very small vector that is not mutually local with R (along qq) can be obtained by acting 
with a finite transformation exp(T^) on R* (along p°). 

In the following sections, we will consider the action of general T-dualities, as we find it 
convenient to describe multi-centre black hole solutions in terms of two auxiliary very small 
vectors R and R* that arise naturally from the equations of motion, as mentioned below 
( |2.6| ). Therefore, we will always consider a T-duality as given explicitly by an exponential 
as in ( p. 48 ), where the explicit action of each order is given by ( 2.35| )-( 2.43D above, rather 



than the equivalent matrix similar to the one in ( |2.3D that has to be defined explicitly. 
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This concludes our discussion of T-dualities in the real basis. In the next section, we 
consider the same transformations in the complex basis, for later use. The reader interested 
in constructing solutions can however safely skip this technical discussion. 

2.3 T-dualities in the complex basis 

In this section we discuss the realisation of T-dualities in the complex basis defined by the 
central charge Z = Z{r) and its Kahler derivative Za = Z(j(r). The discussion here is 
parallel to the one of the previous section, in the real basis, and is complementary to it. 
However, the construction of the T-duality generators in the complex basis will be necessary 
to solve the first order equations describing black hole composites in the following. 

For this purpose, we define the very small vector R as in ( 2.17[ )-( p.l8| ), while R* is 
defined from R using an arbitrary phase, e*° 7^ A/'[J7], as 

Z(i?*) = e3^°/2AA[0]i/2 , Za{R*) = e'''/^M[^f/^Vta , (2.49) 

where the choice of the phases is done for later convenience. It will also be useful to define 
the complex function Y 

_ 2 _ e*"/W[n]V2 _ 21/ 

^ = 1 _ e— AAP - "' Im (e-/Wpi/2) " ^ + '^ ^' ^^.50) 

that has unit real part and the specific parametrisation of the imaginary part will become 



meaningful in the following. The vector (2.49) is by construction mutually nonlocal with 
R because 

{R\R) = (21m (e^^/^^pi/^'jy , (2.51) 



and defines a natural magnetic dual to R. 

We will first determine the T-dualities T+ that leave R invariant. We note that Vta is 
by construction ( p. 181 ) invariant with respect to a subgroup K^ C i^4. In order to describe 
the action of 04 in ( |2.12| ) on this vector, we will parametrize the remaining n^ — 1 generators 
of €4, which describe the coset component €4 €5, in terms of a vector Q". Requiring the 
matrix 

-G\{Q) = CaceC^'^^^dQ" - ^aQ^'Zb , (2.52) 

to be anti-Hermitian and to preserve A/'[Z] fixes the relative coefficients and implies the 
constraints 

QaQ" = , Q" = Minjc^^^nhQc . (2.53) 



Similarly, we parametrize u(l) in ( 2.12 ) by 7 and C"'"' by a complex vector P"^, such that 



the final result we find that the action of Q4 on a general vector reads 
6Z =P^Za + 3i-fZ , 

5Za = PaZ + CabcP^Z' + ilZa + CaceC^'^^^dQ" ^6 " ^aQ^^fe ■ (2.54) 
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In order to describe T-dualities, we must impose that these transformations leave R invari- 
ant, which can be shown to hold if ^ 

p" = -ij^[n]n''-j\f[n]Q^ . (2.55) 

One can now verify that the resulting transformations 

6Z = T+^qZ =3i-fZ-i-fM[n] Cl^Za - M[n] Q^Za 

5Za = T+gZa =i^Za + il M[0] ViaZ-i^ M[n] Cabc^^ Z' - M[n] Q^Z 

- Af[n] CabcQ^Z" + CaceC^'^'^dQ' Z^ - ^aQ^ Z^, , (2.56) 

leave R invariant and commute with each other. Furthermore, one shows that T^g is 
nilpotent of order four, as is clear from the example (^^) where terms at most cubic in the 
parameters A;* appear. Indeed, T^ q can be identified with the corresponding generators in 
(2.35), which act as raising operators on the decomposition ( 2.29| ) and the vector R is the 



highest weight vector, to which we assign weight 3. 

At this point it is important to appreciate the fact that, while we used the complex 
scalar dependent basis to define T-dualities, the following relations hold 

tIqZ{t) = z(T+r) , T^QZaiv) = z„(T+r) , (2.57) 

where T^ denote the representation of these generators in the real basis, parametrised in 
terms of a grade —1 vector A;, as in ( |2.34 )-( 2.35 ). It follows that the parameters 7 and Q" 



depend on the n^ constant parameters k and the scalar fields. 

As alluded to above, the second very small vector (|2.49|) plays a role dual to that of i?, 
as one can check that R* is never a zero mode of the T-duality operator defined in ( 2.5(]| ), 
and in particular 

Z{T+R*) = Z{k) = -3y2AA[0]7, 
Za{TlR*) = Z„(fc) = (3 |y p - 2 y) [7, 7 - 4 i y Q, , (2.58) 

where we also used ( p.37| ) to given the explicit relation of the vector k to the parameters 7 
and Qa- In addition, one can verify that {T't q^R* oc R, as in ( p. 41 ), so that the vector R* 



can be identified with the lowest weight vector of the operators T^n, with assigned weight 
-3. 

As expected from the example in ( p. 3]) , the parameters fc* must be a rank-three vector 
for a general T-duality. Indeed, the vector defined by ( |2.58| ) is of rank three, as can be 
shown by computing its quartic invariant. Furthermore, one can verify that this vector 
satisfies the reality constraint 

Z"" -Minp^'Z = e-'''{d'^^VlbZc + Q.^'iZ -N[n]Q!'Zb)) , (2.59) 



'Note that R is also invariant if 

7 = Q" = , P"n„ = 0, , Pa = -Mp] Caba^^'P" 



These generators, along with the K^ subgroup of K4, and the generators described by (2.55), account for 
the full Gs txR"" subgroup of Ga leaving invariant a given very small vector tell]. 
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introduced in 1 13] in the study of single centre solutions. This defines a Lagrangian subspace 



that includes the small vector R, whereas one also verifies that 

Z{T+R*) = e~2i°Z(T+i?*) , 
n^'ZaiT+R*) = {2e~'''+J^[n])Z{T+R*) , (2.60) 

which implies that {T+R*,R*) = 0. 

Making use of the above structure based on the original very small vector R, one can 
proceed to define similar structures for the dual very small vector R* , in exactly the same 
way. This seems redundant at first sight, since one can always identify the small vector 
invariant under the T-dualities with R, as above. However, this is more natural in view of 
the discussion in the real basis in the previous section, as well as for the applications we 
are interested in, where both vectors appear simultaneously. It is therefore useful to have 
a dual description in terms of R* throughout. 

The T-dualities leaving R* invariant can then be shown to be defined as in ( |2.54 ) with 
parameters given by 

P" = -e'"{ijn'' + Q'') , (2.61) 



where the Q"^ are again constrained by ( |2.53| ), so that the resulting transformations read 

5Z = T-qZ =3i-fZ-i-fe"^ n^Za - e"^ Q'^Za 
6Za = T-QZa = i 7 e-*" naZ-i^ e^" Cabc^'Z'' + i7Za- e-'" QaZ 

- e*" CabcQ'Z^ + CaceC^'^'^dQ" Z^ - ^aQ^' Z^ . (2.62) 

As expected, R is never invariant under ( |2.62 ), with transformation rule 



Z(T,i?) = Z(A;) = -6i — 7 
Za{T-^R) = Za{k) = -2i |y |-2 (y - 2Y)^a 7 + 4 \Y\-^ Qa , (2.63) 

where we show the relation of the vector k in (p. 37) to the parameters in the complex 



basis. Again, one can compute that [T p)^R oc R*, consistent with ( |2.43| ). It follows 
that T~ ^ are lowering operators with R* and R as their lowest and highest weight vectors 
respectively. Furthermore, one can define a Lagrangian subspace that includes R* and 
T~R, through the constraint dual to ( p. 59 ), as 



Z" - e'° n^Z = Mp\c''^''VtbZc + n^ie-'^Z - Q^Zb) . (2.64) 

In addition, the vector ( |2.63| ) satisfies the relations 

Z{T-R) = e-''^M[(l\Z{T-R) , 
^"ZaiT-R) = {2M[^] + e-''^)Z{T-R) , (2.65) 

which imply that {T-R, R) = 0. 

One can obtain the action of the relevant generator hx defined in ( 2.24| ) in the complex 



basis, as the commutator of a T-duality leaving R invariant and a T-duality leaving R* 
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invariant. In general, such a commutator will also give rise to an element of the grade zero 
component 05, as in ( |2.33| ). Choosing the parameters 7 and Q" for the two transformations 
to be identical, one obtains 

[T+,T-] = G(7,g'^), (2.66) 

where 6(7, Q"") is a generator of gl^ 95 bilinear in 7 and Q'^. The gl^ component, which 
is to be identified with the operator hj^, corresponds to the transformation of parameter 
27^ + 3(5"Qa, whereas the 95 transformation is parametrised by 7(5" and Q°'Qh— „ _i {^'i)~ 
■^U,'^VLi))Q^Qc. To project to the gl^ component, one can simply identify the terms quadratic 
in Q" as 

Q'^Qb - ^— {K - l^'^^b) Q'Qc , (2.67) 

so as to cancel all the terms in 95. Equivalently, this identification can be understood as 
an average obtained by acting on the parameter Q'^ with the K^ C K4 subgroup leaving 
Qa invariant, and integrating out the result over K^. By definition, none of the generators 
of gs are K^ singlets, and the resulting expression is necessarily proportional to the gt^ 
generator h^^. Because of the reality constraint (2.53) on the Q"", this average furthermore 
implies 

Q'^Qf' ^^Q^Q<i(mQjc''^c^c - ^n^n^) , (2.68) 

Uy — 1 V ^ ) 

and therefore 

CahcQ^'Q" = ^QbQ' ^f[n]^a , (2.69) 



where we used ( 2.11|) to show that c°-'^'^Cbcd = ^''s^ ^b- ^^ practice, ( ^.69 ) is the only 



constraint one needs to use when computing the commutator in (|2.66| ). 
After imposing the relations above, one finds 

hrZ = ^1^1^(3 (e^"AA[0] - e-'^'Nin]) Z + 2 (AA[f7] - e^") Cl^Zaj , 

hrZa = ^|y|2((e^"AAp - e-'''M[n]) Za + 2 [Mp] - e"*") Zl^^ 



+2 (AA[f7] - e*") Cabcf^^Z^ . (2.70) 



One can now check that (2.44) is indeed satisfied in the complex basis. Similarly, the 
vector space of charges splits into 4 subspaces of eigenvalue {—3, —1, 1, 3} with respect to 
this linear operator. 

By construction, R and R* are the unique vectors of eingenvalue 3 and —3 respectively, 
up to an overall rescaling. The remaining 2n^ directions can then be simply identified with 
the parameters of the two T-dualities T^, as given above. It follows that the eigenspace 



of eigenvalue —1 is spanned by the n^ vectors that satisfy the constraint ( 2.5E ) and are 
mutually local^ with R* , as (i?*,r*'^*) = 0, which is equivalent to 

z(r(-i') = e-2*°z(r(-^)) , 
n''Za{T'~-'^) = (2e-'"+A/'[fi])z(r(-^)) , 

e^^^bZc{T^-^^) = e'^'Z" + e-^"AA[f7]n'^Z(r(-i') . (2.71) 



This is has to be imposed since 71 itself is the trivial solution of (p.59|) 
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Similarly, the eigenspace of eigenvalue 1 is spanned by the n^ vectors that satisfy the dual 
constraint ( 2.64| ) and (i?, F^^') = 0, which lead to 



c'''"'nbZc{r'-''>) = M[n]z^' + n'^ziv'^'^) . (2.72) 

These equations can be identified from ( |2.30| ) in the real basis as ( p. 31 ). 



3. The c*-map, nilpotent orbits and first order systems 

In this section we consider the first order systems describing multi-centre non-BPS black 
holes in A/" = 2 supergravity coupled to n„ vector multiplets labelled by an index i = 
1, . . . ,n„. We refer to the appendix for a short overview of our conventions on 7\A = 2 
supergravity, which coincide with the ones in |13|, to which we refer for further details. 



The systems of black holes we are interested in can be constructed systematically in the 
special case when the special Kahler manifold, A^4, parametrised by the vector multiplet 
scalars, t*, is symmetric. Moreover, we exclusively consider stationary solutions, i.e. we 
always assume a timelike isometry. In this case, one can consider a timelike dimensional 



reduction to three dimensions and dualise all vector fields to scalars |32], to obtain an ef- 
fective euclidean sigma model describing stationary black hole backgrounds. The resulting 
equations of motion are still rather complicated, so that it is common to consider special 
linear systems that solve the full equations of motion, but do not provide a full list of pos- 
sible solutions. There are two such systems known, namely the composite non-BPS system 
1 11] and the almost-BPS system [14|, which together account for a representative majority 
of the explicitly known multi-centre solutions featuring a flat three-dimensional base space. 
The purpose of this section is to define these systems in terms of four-dimensional, mani- 
festly duality covariant quantities, aiming for a clear description of their general structure. 
To this end, we make use of the duality symmetries of the three-dimensional theory 



resulting from the dimensional reduction, through the formalism developed in [33, [I 



We therefore first describe the basics of this effective theory in section 3.1, followed by a 
discussion of the method of nilpotent orbits in section |3.2| . We then present the derivation 
of the first order flow equations for the composite non-BPS system and the almost-BPS 



system in sections |3.3| and 3.4 respectively. Note that, while the derivation of the almost- 



BPS system has not appeared before, our section 3.3 is essentially a review of the derivation 



of the same system in [13|, which we include for completeness. 



3.1 The three-dimensional non-linear sigma model 

In order to describe stationary asymptotically flat extremal black holes, we introduce the 
standard Ansatz for the metric 

ds2^-e2^((it + a;)2 + e-2^dx-dx, (3.1) 

in terms of a scale function U{x) and the Kaluza-Klein one-form uj{x) (with spatial com- 
ponents only), which are both required to asymptote to zero at spatial infinity. Here and 
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henceforth, aU quantities are independent of time, so that ah scalars and forms are de- 
fined on the flat three-dimensional base. The n^ + 1 gauge fields of the theory F = dA 
for / = 0, . . . n^ include the graviphoton and the vector multiplet gauge fields. Together 
with their magnetic duals, these can be arranged in a symplectic vector, J-", as in ( [A.2| ), 
which transforms linearly under electric/magnetic duality. For a background as in (p.lj), 
the appropriate decomposition of the gauge fields takes the form 

2A = C{dt + (jj) + w (3.2) 

and accordingly for the field strengths 

2T = dC{dt + uj)+F, F = Cduj + dw, (3.3) 

where we defined the gauge field scalars C, arising as the time component of the gauge fields, 
and the one-forms w describing the charges. Here, F is defined as the spatial component 
of the field strength, which is not closed but satisfies 

dF = e^^ • dcj A JF , (3.4) 

according to ( [A. 15 ), which can be written as 

d^ = e2^J*F. (3.5) 

Note that this first order equation determines the C, in terms of the vector fields w and the 
scalars. 

Upon dimensional reduction over the time direction, the set of moduli i* parametrising 
the coset space (|2.7D are extended to include the scaling factor U and the scalar dual to the 
angular momentum w in ( |3.lD , as well as the fields C, in ( |3.3D , which altogether parametrize 
the para-quaternionic symmetric space ^^ 

Ms = Gs/{SL{2) X G4) . (3.6) 



This defines the so-called c*-map, which can be related to the standard c-map |3^ by 
analytic continuation. The three-dimensional symmetry group Lie algebra 03 decomposes 
as 

03 - 1(-2) e [(-^' e (0[, e 04)^ e \'l' e 1'^' , (3.7) 

where the weights refer to the eigenvalues under the adjoint action of the q\i generator. 
The grade one generators in v^^ are associated to the gauge invariance with respect to a 
constant shift of the scalars C,, and accordingly the grade two generator corresponds to the 
shift of integration constant defining the scalar dual to w. 

As discussed in more detail in p3[] , the equations of motion for the scalar fields parame- 
trising the symmetric space M.^ are expressed in terms of the corresponding Maurer-Cartan 
form 

v-^dv = P + B. (3.8) 



^"Here para-quaternionic refers to the property that the holonomy group of Ms SL{2) x G4 C 51/(2) x 
Sp{2n^ +2,R). 
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Here, v £ G3 is a coset representative describing the scalar fields, while P is the coset 
component of the Maurer-Cartan form, defining the scalar momenta. Similarly, B is the 
sb © fl4 component defining the pulled back spin connection. In components, the scalar 
momenta are defined as 

w = -dU - ^-e^^ * dw , S" = -e^df , Z = e^Z(*F) , Za = e^ Za{-kF) , (3.9) 



where we introduce some shorthand notations, based on the central charges in ( A.ll| ), that 
will be used for the remainder of the section. 

At this stage it is important to introduce some properties of the 03 algebra. The 
components of an element of the Lie algebra 93 in the coset component 93 (s[2 © 54), as 
given in (|3.9| ) above, correspond to U{1) x U{1) x i^4 irreducible representations as the 
two complex parameters w and Z and the two complex vectors Z"',Ti°' transform in the 



C""" representation of K^^, according to the decomposition ( 2.12 ), the same as the scalar 
field momenta e^dt* in four dimensions. As one can check explicitly, the quadratic trace 
invariant which defines the SL{2) x G4 invariant norm 

|w|2_ iZp-Z^Z'^ + S^Sa , (3.10) 



is equivalent to the effective Lagrangian in the background ( |3.lD . 
The 5(2 algebra is realised on these components as 

^w = ipw + AZ, 5Z = -ipZ + Xw, SZ"" = ipZ'' + \T.^ , ^S'^ = -ipS" + AZ'* , (3.11) 

where p and A are a real and a complex parameter respectively, parametrising the s[2 group. 
Similarly, the action of 94 can be written as using the parameters defined below ( p.l2|) , as 

(5w = KZ" + 3i7w , 5S" = A"Z + c^^^AfeZ^ + G"feS* + ^7^" , 
5Z = Kj:" + 3«7Z , dZ"" = A"w + d'^^Kb^c + C'bZ^ + i-^Z"" . 

Note that the 94 action defined by these equations corresponds to the divisor group in (^), 
rather than the original 94 as given in four dimensions in ( p. 13 ), the two being related by 
a conjugation in 93. 

Finally, we give the components of the SI2 ©94 component of the Maurer-Cartan form, 
B, along s[2 

p[B) = -\e^^ i. duj - \Q , \{B) = e^Z{i.F), (3.13) 

and along 94 

^{B) = -\e''^^duj + \Q, Aa{B) = e^ Zai^F) , (3.14) 

CbiB) = e^djei dP - ajbdie^^df - -6^Q , (3.15) 

where G"'b{B) defines the €4 valued traceless^^ component of the pulled back spin connection 
on M-4 and Q is the pulled back Kahler connection ([A.8| ). 



^To prove that G"i,(i3) is indeed traceless, one can use (2.1C) and (^.ll[) to show that c'^'"''Cbcd = ^^^^4-^5^. 
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3.2 Nilpotent orbits 

The basic observation for constructing black hole solutions in four dimensions using the 
three-dimensional Euclidean theory describing stationary solutions is that regular station- 
ary solutions of A/" = 2 supergravity with a flat three-dimensional base metric can be 
described by a three-dimensional momentum P that is nilpotent as a Lie algebra element. 
This implies in particular that P can be written in terms of the basis element e^ of a nilpo- 
tent subalgebra of 03. Such a subalgebra is always associated to a semi-simple element ^^ 
h of s [2 © 04 such that 

hea := [h,eQ] =p„eQ, , 1 < Po < n , (3.16) 

where n defines the maximal possible eigenvalue of adh in 03. This implies for instance the 
equation 

n 

n(h-i)P = 0, (3.17) 



which defines a first order constraint on the components of P. In order for ( p. 17 ) to be 



consistent with the equations of motion and the Bianchi identity, the covariant derivative 
of the generator h must satisfy 

n n i— 1 n n i—1 

Y, n (h-i)dBhJ](h-fc)AP = 0, Yl n (h-J>BhJ](h-A:)*P = 0. (3.18) 

i=l j=i+l k=l i=l j=i+l fc=l 

These equations are satisfied if dsh also lies in the nilpotent algebra defined by h, or 
equivalently if one imposes that 

n 

JJ(adh - i)dBh = . (3.19) 

1=1 

Note that this is not necessarily the most general solution for the generator h, as it might 
be possible to construct special solutions, for which e.g. there are preferred spacetime 
directions described by a nonzero derivative. Considering this type of solutions can be 
understood as the natural generalisation of the BPS black hole solutions in [^, ^ to include 
supersymmetric string solutions. These are not accounted for by the standard black hole 
ansatz, but are are allowed by the requirement of preserved supersymmetry. Similarly, we 
will only consider the generic situation in ( 3.19| ), which should be satisfied for all composite 



black hole solutions with a flat three-dimensional base. 

Without loss of generality, one can always choose the generators h G 5(2 04 such 
that only its components A and Aq do not vanish, as these are the only generators with a 
positive Cartan norm^'^. Equation (3.19) can then be viewed as first order equations for 
these auxiliary components, which can be solved to determine their evolution in space in 
terms of the physical fields. With this information at hand, equation ( p. 17 ) then defines 



^^Semi-simple means that it is in the conjugation class of an element of the Cartan subalgebra. 
^^This is similar to the canonical choice of d/dt as a timelike vector in Minkowski spacetime, i.e. there 
are more general timelike vectors but they can be obtained from the canonical one by symmetries. 
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first order equations for tlie pliysical fields, wliich contain the auxiliary components A and 
Aq and determine w and S" of ( p.S| ) in terms of e Z{-kF) and e Za{*F), plus some possible 
constraints on the latter if the dimension of the coset component of the nilpotent algebra 
defined by h is strictly less than 2nj, + 2. 

Given the definitions in the previous section, it is possible to make ( 3.17 )-( ^?l9D explicit, 



in order to determine the auxiliary components and the first order flow directly. A simple 
example is given by the BPS system, which is characterised by an element of h* Gsb, i-G. 
A(j(h*) = 0. It follows that its only nonvanishing component is A(h^,) = e*°, where this 
phase is identified with the phase a that defines the covariantly constant spinors as in Q. 
The action of this generator follows from ( ^.11 ) as 



h,w = e-*"Z, h,Z = e^"w, h,Z'' = e-^"S", h,S'^ = e^"Z", (3.20) 

while the relevant eigenvalue equation ( p.l7| ) is simply 

KP = P- (3.21) 

The last relation clearly imposes a linear relation between the derivatives of the four- 
dimensional scalars and the gauge fields upon using ( 3.20| ). The phase e*" is determined 
by ( |3.19| ) , which can be shown to reduce to 

da + Q + ]-e^^ • dw + 2 Im(e-^°Z) = da + Q - ]-e^^ * dw = , (3.22) 



where in the second equality we used the first of ( ^.20D and ( ^.21 ). These equations are 



easily seen to be equivalent to the BPS system of ||7[ ^ and we refer to ||3^ for a more 



detailed analysis of the eigenvalue equation ( 3.16 ) for this generator, leading to the system 



of equations describing multi-centre BPS black holes. Similarly, non-BPS solutions with 
vanishing central charge at the horizons are described in a similar fashion, with A = and 
a normalised rank one A^, (i.e. c'*'"^Af,Ac = and AqA" = 1) p^ . 

Other examples of such constructions using nilpotent orbits have been used to obtain 
the systems describing respectively single centre non-BPS black holes with a non-vanishing 



central charge at the horizon and the composite non-BPS multi-centre system |13]. In 
the next two sections we discuss in some detail the construction of the two multi-centre 
non-BPS systems, the composite non-BPS and the almost-BPS, which form the basis of 
this paper. 

3.3 Composite non-BPS flows 

The composite non-BPS system |11, ^, 13 1, describes configurations of interacting non- 



BPS centres and corresponds to an eigenvalue equation as in ( |3.16| ), where the relevant 
element he belongs to 04. As the scalar momentum P must be of positive grade with 
respect to he (given that pa > 0), we consider the relevant decomposition of the coset 
component, parametrising 2 ® ]]^2ni,+2 through ( |3.9| ), in terms of grades with respect to an 
element he € 34. The relevant graded decomposition is in this case 

04 = (R"")(-^) e (sii e 05)*"' e (R"")<'' , (3.23) 
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for 04 itself and 



2 (gl R2n„+2 ^ 2(-3) 



n„\{-i) 



R"") 



iE 



n„N(i) 



)(3) 



(3.24) 



for the coset component, i.e. P E (2 (8) R"")'^^ © 2'^'. We can choose he to be Hermitian 
(i.e. to he in 94 (u(l) ^4)), so that it is reahsed for 



7(hc)=0, Aa(hc 



n. 



G\{ha)=0, 



(3.25) 



where Q,a describes a very small vector, as in ( p.l8| ). Equivalently, ^a is in the [/(I) x K4 
orbit of the Jordan algebra identity. 

More explicitly, one finds the following action on the coset component 



hcw = 0„Z" hcZ = 0„S" hcZ" = S7"w + c'^'^Qb^c hcS'^ = Cl^Z + c'^'^nhZc , (3.26) 

which is to be identified with the elements of the coset component, through ( |3.16| ) for 
Pa = 1 and Pa = 3, according to (p. 24). Considering a general linear combination of the 



grade one and three solutions as in |13], one can express w and S" in terms of Z and Za as 

^abcri'y I rtafry \ /Toi nt"7_ \ (3 27) 



w = -{VLaZ^ - M[n]Z) , S" = c^^^VLbZc + -rj"(Z - M[n]VL''Zb) , 



which are the explicit first order relations for the scalar momenta, upon using the definitions 
in dp). 

These flow equations contain the auxiliary components N\yi\ and ila of he, which 



satisfy ( 2.181 ) and define a very small vector R of unit mass through 

Z{R)=M[n] Za{R) = VLa. 



(3.28) 



The flow equations for these auxiliary flelds are given by ( 3.19| ), which in this system 
reduces to 

[hc,dBhc] = 2dBhc , (3.29) 

in view the fact that d^hc lies in 34 and is inert under 95 by definition, and the decom- 
position in ( p. 23 ). Using the explicit form of B in ( p.l3 )-(3.15) and the first order flow 
(3.27), one computes the components of d^h as 



Ka{dBl^) = Za{dR) + MMp]Vtbe\df]Za{R) -(Za+ M[n]Cabc^' Z' - fl, Vi^Z 



G\{dB\y) = c'^^'cMei^'^Zc - nj'') + QbZ" - n'^Zb - ^6^1m[Q'Z, 



(3.30) 



where we explicitly separated the terms depending on the derivative of the vector R. It is 
now straightforward (though cumbersome) to evaluate ( 3.2g| ), imposing the above relations 
and 

c''^^Cbde^'^{J^['^]Ccfg^^ Z>^ - ^c^fZ^) - fl'' {f^[n]cbcd^^ Z'^ - nb^cZ"") 

= -d'^'Cbde^cZ'^ + ^bZ" , (3.31) 
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which can be shown using ( ^.11 ). The result one finds is that ( |3.29| ) is identically satisfied 



for ( 3.3C| ), if and only if the combination of the first two terms in Aa(dBh) vanishes. This 



condition implies |13] that there exists a constant symplectic vector R, since 



dR = -Re[M[n]nhe''idf]R , ^ R= ^. , (3.32) 

\Z{R)\ 

where the second equation follows from the first by use of standard special geometry iden- 
tities. It follows that the generator he is in this case determined by a constant very small 
projective vector R and the scalar fields are such that 

A.(h.).||. (3.33) 



One can now return to ( 3.27| ), which becomes a first order flow equation for the scalars dU, 



-kdw and dt* in terms of the gauge fields and the constant vector R. Further details and 
the characterisation of solutions to these equations are given in section ^ 

A special case of this system, namely the single centre subclass, was discussed in 
detail in ||l^. Indeed, one expects a system describing multi-centre black hole solutions to 
contain a consistent subsystem describing single centre black holes. In the case at hand, 
one can restrict the above equations ( p.27| ) to this subsystem by imposing that the scalar 
momentum is also of positive grade with respect to the generator i(hc -|- h*) where he 
and h* are the generators that define the composite system through (|3.26D and the BPS 
system through ( 3.20| ) respectively. The result one obtains is that the single centre non- 



BPS momenta satisfy (|3.27 ), for Z and Za constrained to satisfy exactly the same phase 



dependent projection ( 2.59| ) above, 

Z" -M[^]Q.''Z = e-''^{c''^''VLi,Zc + 9.''{Z -M[^]0!'Zb)) , (3.34) 

which represents a constraint on the physical degrees of freedom that is necessary to reduce 
to single centre solutions. One can interpret this constraint in terms of the decomposition 
in ( 2.29 ), upon identifying the very small vector R in (|3.33D with the vector R in section ^^21 . 



It follows that ( |3.34 ) implies that only the grade (—1) and (-1-3) components of the gauge 



fields (and thus charges) are allowed. We refrain from discussing the single centre system, 
as we will deal with the multi-centre case in what follows. However, we will sometimes 
make use of structures present already in the single centre case, referring to |13] for an in 
depth discussion. 

3.4 Almost-BPS flows 

We now turn to the derivation of the first order fiow equations for the almost-BPS system 



1 14] of multi-centre black hole solutions. Our treatment is purely in terms of four /three 
dimensional quantities, but we still refer to this system as almost-BPS, as we will find 
that it contains the original almost-BPS equations, derived from the BPS conditions of 
five dimensional supergravity. However, the five dimensional system is only defined in a 
fixed frame and does not allow for generic charges at the centres, while we will present a 
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manifestly duality covariant form of these equations, as we did for the composite non-BPS 
system in the previous section. In this sense, the almost-BPS system below is larger, as 
it contains the one in five dimensions, as well as all the charge configurations that can be 
obtained from it by dualities. Our covariant system also represents the general form of the 
almost-BPS system of [^], where a frame was chosen to match with the five dimensional 
results. Closely related first order flow equations have recently been derived using similar 
methods in |^] for the STU model. ^^ In this case Q.a reduces to three phases, which 
altogether with a define equivalently the first order equations. Nevertheless, the authors 



in 1 21] use the second order equations of motion to determine the evolution of these phases, 
whereas we consider the first order equations ( p. 19 ), which imply the second order ones by 
construction. 

For the almost-BPS system, the eigenvalue equation ( |3.16| ) is defined by the generator 
Ha = he + 2h*, where we use the same generators as defined in (3.26) and (p.20|) above. 



The corresponding graded decomposition of s[2 ® 04 is then given by 

512 ® 04 = 1'-'' e (]R"")(-^' (0ii egti ®05)'°' e OR^'" ® i<'' , (3.35) 

which is exactly the same as the one in ( |3.23| ) with respect to the 04 components. The 
additional notation R"" and R"" refers to two conjugate representations of the real group 
05. Similarly, the coset component is decomposed as 

2 ® R2n„+2 ^ -,_(-5) ^ (J^ri^y-S) (1 ]R".)(-i) (1 R^")*^' (R"-)'^' ® l'^' , (3.36) 

so that P G (1 © R"")<^' © (R""')<^' © 1<^', i.e. it lies in the positive grade component. 

Due to the presence of the generator h*, which by itself defines the BPS system, ^^ the 
first order system reduces to the BPS system (p. 21) for most of the components of the field 



strength F, except for the ones that define the solution 

Z = iet°A/'P^u, Z„ = ze5"AA[n]5 0aU, (3.37) 



for which the expressions of w and S" take the opposite sign, where u is a real function. 
The projection to this component is obtained as 

VZa = J e^Np]^ J^almfe-'^AAp^Z + e'^f M[n]^Tt Z^] , (3.38) 

for which one easily computes that "P^ = V. One can then verify that the first order system 



"i.e. for Mi ^ SU{1, l)/U{l) x 5(7(1, 1)/U{1) x S'f/(1, l)/U{l). 

^""Note that while the same generator appears in the reduction of the composite non-BPS system to the 



single centre class, we treat Ha as the defining object here, instead of imposing a constraint similar to (2.59) 
on an existing system. 
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obtained by changing the sign of the components above, given by 



w 



'{Z - 2VZ) 



l^e-^-Z + \MmZ + \^aZ'^ -\e- 



W[f7]0"z, , 



^a ^ e'^^iZ" -2VZ'') 



e^^Z"" + -^^ {Z + M[Q]d' Zb - e*" 



M[n]Z - e'^'UbZ'') , 



(3.39) 



solves ( 3.16 ), as expected. Using the definitions in (3^), these equations become a first 
order flow for the four dimensional scalars in terms of the gauge field strengths. 

However, the auxiliary fields e~*" and f^a parametrising hA are still arbitrary, we 
therefore need to impose ( ^.19 ) to obtain their dependence on the physical fields. First, 
the condition on dsh* to be of positive grade is the same as the one for the BPS system 



in ( 3.22 ) and leads to the condition, 



da + Q + -e^^ i.doj + 2 9[e" 



'Z] = 0, 



(3.40) 



where the last term is now given in terms of scalars by ( p.39| ) and is therefore more compli- 
cated than in ( ^.22 ). Analysing the components of ( |3.19| ) along 04 to obtain the remaining 
auxiliary components, $7^, one finds that it reduces again to ( [3.29 ) in exactly the same 
way as for the corresponding quantities in the composite non-BPS system. In fact, one 
can easily check that ^{(Ib^a) and G°'i){dB^A) are exactly the same as in ( 3.30| ) since they 
only depend on the algebra ( p. 13 ). In contrast, manipulating the expression for Ka{dB^A) 
is sensitive to the particular form for the scalar momenta in terms of the gauge fields. In 
order to write the resulting equation in a suggestive way, it is useful to draw intuition 
from the analogous vector in the single centre class^^, as well as from known multi-centre 
solutions, which indicate that the magnetic dual of the vector described by Via is a more 
convenient variable. We therefore consider the dual very small vector of mass one, i?*, 
defined as in (|2.4g| ), i.e. we impose 



0, = e'''Z{R*)Za{R*) . 



(3.41) 



With this definition, one can use (3.39) to show that 



K{dB^A) = e^" [Z{dR*)Za{R*) + Z{R*)Za{dR*) + 2Re[e*"flfee^(if ] Z{R*)Za{R*) 



Za + Af[n]cabc^'z''-nanhZ'). 



(3.42) 



As in ( |3.30| ), the second bracket in the last expression satisfies ( 3.29] ) identically, whereas 
the first term vanishes upon imposing that R* is related to a constant vector R* by 



dR* = -Re[Z{R*)ZbiR*)e''idf]R* 



R* 



R* 



\Z{R*) 



(3.43) 



^^Note that the non-BPS single centre system is also contained in the almost-BPS system, but we refrain 
from discussing the details of this reduction. 
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similar to ( ^.32 ) for the composite system. In contrast to ( |3.3[l| ), there are two terms 



containing the derivative of R* in (3.42), so that one could expect more general solutions 
than the one shown above. However, one can verify that any solution for R* other than 
( |3.43 ) is such that it is mutually nonlocal with its derivative dR*, which is not allowed for 
an everywhere very small vector. 

This concludes our presentation of the derivation of the first order flow equations for 
the almost-BPS system. One may now consider solutions to (3.3£), which is a first order 
flow equation for the scalars dU, -kduj and dt^ in terms of the gauge fields and the constant 
vector R* , upon using ( 3.41| ) and ( p. 43 ). In section ^ we discuss the real form of these 
equations, show that they correspond to a linear system and give the characterisation of 
their solutions in terms of local functions. 



4. Composite non-BPS system 



In this section, we present in detail the steps required to characterise solutions of the flow 
equations for the composite non-BPS system in terms of local functions. The starting point 
is the solution of the nilpotency condition ( 3.27| ), written explicitly as a first order flow 
system for the four dimensional scalars and the metric degrees of freedom 



dC/ + -e2^ 



• duo 
efdf 



--e^iQaZ'^i^F) -J\f[n]Z{*F)) 

c"^^Jlfee^Zc(*F) + i J7"e^ {Zi*F) - M[n]n^ Zb{*F)) , (4.1) 



where F is the spatial component of the field strengths defined in ( |3.3|) . The vector JIq is 
related to the constant very small vector, R, through ( ^.28 ) and ( p. 32 ) above. For later 
reference, we give the inverse relations for the field strengths 



e^Z(*F) = -M[n](dU 



2 



1 



nbe°df 



'Za{*F) 



-Cabc^^eldf + -Mp\VLa^be\df - -VtaidU - -e^^ * duj 



(4.2) 



In order to solve this system, we first construct the electromagnetic potentials and use 
the resulting structure to simplify the equations in section 4.1. In section |4.2| we exhibit 
the relevance of the T-dualities introduced in section ^. We then proceed to rewrite these 
equations as a linear system of differential equations and discuss its integration in terms of 
local functions in sections ^^ and 4.4 respectively. The reader interested in applications 
can find a summary of the final form of the system in section [4.5|. 



4.1 The electromagnetic potentials 



In order to solve the system (4.1), we construct the gauge field momenta ( |3.5| ) using the 
derivative of R given by ( 3.32] ), to obtain 



dC = 2e^^K[Z(*F)V + Z*(*F)AV] 



d {e^^M[n]V - n'DiV]) + ie^ ^MinjQidf] R 



\e^^i^dujR, (4.3) 
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where we made extensive use of the special geometry identities in section ^. The first term 
is manifestly a total derivative, whereas the others are along the very small vector R and 
must therefore combine into the derivative of a single function. This requirement, along 
with (3.32), leads to the condition 



e^ ^[Np] n,df] + ie^^ i.doj = M e"^ ^{M{VL\^idt\ - d{M e^^) , (4.4) 

where M is an arbitrary function, so that the gauge field momenta take the form 

C = e^^\M{VL\V - 0^ AV] - )-e'^MR , (4.5) 

with the corresponding central charges given by 

Z{Q = ^ e^ (1 + i e2^M)AA[f7] , Z,(C) = ^ e^ (1 + i e^^M)^^, , (4.6) 

for later reference. 

The structure of (|4.5| ) can be used to show that one vector is always trivial, simplifying 
the system. To see that, we compute 

(^,0 = -2e^|^(i?)|, (4.7) 

whereas taking the imaginary part of (^]^) one gets 

e2^*da; = -^- {R,i.F) . (4.8) 

2\Z{R)\ 

Finally, using (^) , one finds that 

{R,dw) = Q, (4.9) 

which implies that one vector field is always absent. In terms of the graded decomposition 
in section 0, the vanishing component is along the very small vector dual to R, given in 



(2.49), as will be shown shortly. 

One can now combine (^^) and ( |4.9| ) to disentangle the term proportional to -kduj in 



the definition of the scalar fiow equation, so that ( |4.1| ) becomes 

- eldf = d'^'^ni^e^ Zc{i^dw) + \Q.''e^ {Z{i.dw) - M[^]Q^ Zi,{irdw)) 

+\e^^{-M + ie-'^^)N[n]Q.''i.duj. (4.10) 

Applying the same procedure on ( |4.2D , one obtains the inverse relations of ( [4.10| ) for the 
central charges Z{dw) and Za{dw). The charge vectors dw can then be straightforwardly 
constructed with the result 

i.dw = 2'^[-Z{dw)V + Z''{dw)DaV] 

= e-^Q[-{dU + fe^^ • duj)n+ + {iY e^^ * do; - M[n]n'dti)n^ 

-2 c'^^^nbeldPDaV + 2J\f[n]n'di0DjV] , (4.11) 



31 



where we used the shorthands 



7^± = ± M[n]v + n'DiV , 

Y = {l + ie^^M). 



(4.12) 
(4.13) 



Note that the above results are in direct correspondence with the ones in [13|, where the 
single centre system was treated. 

4.2 Connection to T-dualities 

Given the structure above, it is useful to define a second distinguished very small vector 
which is mutually nonlocal with i?, as in (|2.49|)-(p.50) through 



R* 



Z(i?)|"^Re Y'^ N[^]V + \Y\^YQ}Diy 



(4.14) 



Here, we identified the function Y in ( [4. 131) with the one in ( 2.50 ) and we included an overall 
rescaling. One can check that (i?, R*) = A and that the central charges of R* above satisfy 



(2.18). This new vector is defined in exactly the same way as the second constant vector 



used in |13] in the single centre case, but is not constant in the full composite non-BPS 
system. 

As explained in (p.5|) = (p^), a non-constant very small vector with a constant non-zero 
symplectic product with R must be of the form 



R* = exp[T+]i?;5 



(4.15) 



where Rq is a constant very small vector satisfying {R, Rq) = 4 and exp[T+] is a T-duality 
matrix leaving R invariant. It then follows that the derivative of R* takes the form 



(IR* -dT+R* = 0, 



(4.16) 



which is identically closed by the fact that T-dualities are abelian. The above are consistent 
with known solutions [|ll| , [l^ , in which R* takes the form of a T-duality whose parameters 
are harmonic functions, acting on a constant vector along p^. 

Within the composite non-BPS system, one can explicitly compute the components of 
the derivative of R* using ( p^ ) and ( |4.10[ ) , as 



Z{dR*) =T+qZ{R*) = Z{dT+R*) , 
ZaidR*) =T+Za{R*) = ZaidT+R* 



(4.17) 



where we indicated that the components in this basis are given by the variation of R* itself 
under the T-duality transformations ( 2.56[ ), as shown in (2.58). In order to obtain this 
result, one has to identify R as the grade (+3) very small vector, that is invariant under 
all T-dualities T+, which we assume henceforth. The explicit values for the parameters are 
given by 



7+ 



1 21/ 
6 ^ 



{dM - -kdoj) , 

^[YMMNa + YCabc^'^N' 



sta (2Y Af[n]nbN'' + Y M[n]n''Nb^] , (4.18) 
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where we used the combination 



N"" = -e^df + M[n] n^ldU- ^-e^^ * do; ) , (4.19) 



for brevity. Note that in the single center class of [13| these expressions were shown to 
vanish, consistent with the fact that R* is constant if the T-dualities are rigid. 

Given the above, we can directly apply all considerations of section |2|, since the presence 



of the two constant very small vectors R and Rq implies that the grading shown in ( 2.2(]| ) 



is relevant for the integration of the system. It follows that Rq is identified as the grade 
(—3) very small vector. As discussed in ( 2.58D -( ^l59| ), a generic T-duality is parametrised 



by a rank three grade (+1) vector, which we denote by /C, so that ( 4.15| )-(4.16) become 



R* = exp[T+]i?S , dT+^Rl = dJC . (4.20) 

The rank of /C can be verified by checking that the quartic invariant of the vector dR* in 
(4.17) is vanishing, i.e. 



hidR*) = h{dT+R*o) = h{d)C) = , (4.21) 

d^T+d,T+R*o = -^ li{d^T^R*Q, d,T+R*Q, R) = -^ l'^{d^)C, d,)C, R) . (4.22) 

In what follows, we will generally drop the subscript /C on T^ for simplicity, since this is 
the only T-duality appearing throughout. 

Finally, it is worth commenting on the difference between the dual very small vector 
in ( 2.49|) , which was used to define R* in ( 4.14| ), and the constant vector R^ of grade (—3), 



that might seem confusing. It is important to realise that (|2.49| ) simply defines a possible 
dual vector, which is not unique. Since two very small vectors not commuting with R 
are related by exactly a finite T-duality leaving R invariant, one may choose any other 
vector in that orbit. We will fix this ambiguity by defining the function /C to vanish in the 
asymptotic region, or equivalently Rq = R*\r^oo- 

4.3 The linear system 



One can now use the vector R* of ( 4.14 ) in a way similar to the vector R was used in 
section 4.1, to project the flow equations and simplify the system. Indeed, taking the inner 
product with dw one finds 

{R*,dw) = -*(i(e-^|Z(^)|-i|yp) = -^dV, (4.23) 

where we defined the function V. Note that V is not a harmonic function, since R* is 
not a constant vector. Combining this with ( |3.32| ), ( [1.4| ) and ( [4.101 ), we can determine the 
combination Ojdi' in terms of V, M and the metric components as 

^M[n]nidf] = \z{R)\-'^d\z{R)\ = -v-^dv - du + \y\-^ d\Y\^ , 

y2 
J\f[n]ni df = (1 - 2Y) dU - \i e^^ i. duj + —— dV - V'^Y dV . (4.24) 
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We now use all the above information to write the expression ( [4.11 ) for dw in a sugges- 
tive form. To this end, we use the expression for the derivative of R* in ( 4.17] ), combined 



with inspiration drawn from the analogous computation performed in [I3| for the single 
centre class. After a long but straightforward computation we obtain 

-kdw 



Ae-^lm[3iYAf[n]V - i{2Y - Y)n'DiV] 7+ 



+e-^ ^[(2 dU - i e^^ • duj)e-"'V - 2 e-'^'dfDiV - 4 Q^ DaV] . (4.25) 



Here, we used ( 2.50 ) to define a new phase as 






M[n] 



(4.26) 



which will be useful in what follows. One can compute the derivative of e *° using (|3.3 
to first obtain 



Af[n]id-iQ)M[n] - 2 

from which follows the relation 



^e^^ + dw + ^^dy 



1 2(7 



IrP 



2U i/M- 



le 



2(7 



Vd{^) 



da + Q + ie^^ i^doj- e^^ d(f-) V = 0. 



(4.27) 



(4.28) 



We emphasise that these equations are completely analogous to the ones relevant for the 
single centre system, but now involve generically non-harmonic functions M and V. 

Additionally, ( 4.25| ) contains the parameters for the particular T-duality appearing in 
the derivative of R* in ( |4.18 ). These can be rewritten by observing that the action of a 
T-duality on the symplectic section takes the form 

dT+[2e-^Im(e-*"V) ^^ " 



V 



R* 



= 2 e"^ Im [-3 i 7* y e-^° V + i 7, (y - 2 y)ff AV - Qt A V] 
Using this relation and ( [4. 28] ) in ( [4.25| ) leads to 

i^dw = - [d-2dT+] [2Im(e-^-^°V) 



^VR 



M £>* 

■yR 



(4.29) 



(4.30) 



where we also used the fact that R is by definition inert under the T-dualities in question. 
This is the final form of the flow equations in the real basis, where the local T-dualities 



have parameters given by ( 4.18| ) above. 

4.4 Integration and local structure 

Due to the presence of the flat connection for the T-dualities, it is not possible to solve the 



system of equations (4.30) in terms of harmonic functions only. However, the scalar and 
vector fields can be written as 



2e-^Im{e-'"V) - ^V R - y 



M jj* 

dw = {-kd-2 i.dT+)'H, 



-n, 



(4.31) 
(4.32) 
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where the vector of functions % is the solution to the non-harmonic equation 

dirdn-2d*dT+n-2dT+A^dn = 0. (4.33) 

In order to disentangle the derivatives on T-l and T+ and cast this as a Poisson equation, 
we introduce the rescaled vector of functions 

?^o = exp[-T+]^, (4.34) 

in terms of which we find 

2 e-^Im(e-^°V) -^VR-fR* = - exp[T+]Ho , (4.35) 

dw = exp[T+] * dJio - W(exp[T+]) Tio , (4.36) 

while ( 4.33| ) takes the form 

di^dUo-di^ dT+Uo - dT+ A i^dT+Uo = . (4.37) 

It is possible to give a systematic characterisation of the solution for T-Lq and exp[T+] 
using the following crucial observation. From the expression ( 4.36| ) for the vector fields, we 
compute for the derivative of Hq that 

Za{exp[T+] • dUo) - J\fM ^a Z(exp[T+] * dTio) 

= Za{-kdw + exp[T+]dT+?^o) - ^^P] ^a Z{-kdw + exp[T+]dT+?^o) 

^ {Q+a + 2iQa7+)+e'''Na, (4.38) 



ll'l 



where Na is the combination of scalar momenta defined in ( [4.191 ). One can easily show 
that the central charges in ( 4.38 ) satisfy the reality constraint (|2.59|) , i.e. we have 



c'''"=S7bZc(exp[T+] d-Ho) + n"" f Z(exp[T+] dTio) - f^M n^Zb{exp[T+] dUo) ] , (4.39) 



This constraint was found to be crucial to describe single centre solutions in [^], where 
it was analysed in some detail. Using ( p.30|) , one computes that it is equivalent to the 
constraints 

]-l'i{R,R\ex^[T+]d'Ho) = -4exp[T+]d?^o + 3(exp[T+](i?^o,i^*)^ , 

(i2,exp[T+](i?^o) = . (4.40) 

These equations are manifestly duality covariant, and in particular T-duality covariant. 



Therefore, using ([4.15) one finds that the vector dl-Lo satisfies the constraints 



h'i{R, Ro, dUo) = -WHo + 3{dno,R*o)R > (^> '^^o) = . (4.41) 

In a similar fashion, one can show that the vector of functions T-Lq itself satisfies the same 
constraint, using ( |4.31 ), and therefore only half of its components are allowed. 
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Indeed, a vector satisfying this constraint does not contain components of grade (+1) 
and (—3), so that in terms of ( |2.2E ), one finds 



Ho € (R"-)(-^'eE('' , (4.42) 

and therefore describes only n^, + 1 functions instead of the 2 [riy + 1) one would have a 
priori, i.e. it lies on a Lagrangian subspace. One of these functions is clearly the function 
V in ( [4.23 ), describing the grade (+3) component of T-Lq along the direction of R. The 



remaining n^ functions span the grade (—1) vector space, according to the decomposition 
imposed by the T-duality T+, and are undetermined for the moment. 

Using this decomposition in ( [4.37 ), we find the following grade assignments for each 



term 

di^dUo-d* dT+no - dT+ A ^dT+Tio = , (4.43) 

(-l)®(+3) (+1) (+3) 

where we used the fact that T+ is of grade (+2). Now, since the second term in ( [4. 43)) lies 
in a subspace orthogonal to both the other terms, it is clear that this equation decomposes 
in two independent equations, as 

di^dno-dT+A*dT+no = 0, (4.44) 

d*dT+ = 0. (4.45) 

In deriving the second equation we used the fact that applying a T-duality on a generic 
vector as in ( |4.42| ) results in a vector of grade (+1), which cannot vanish for a physical 
solution^^ unless the T-duality matrix is trivial. We now analyse each of the two equations 
in turn. 



First, it is easy to express (4.45) in terms of the vector of T-duality parameters in 
(g^, as 

d*d/C = 0, (4.46) 

so that the parameters fC are identified as a rank three vector of harmonic functions. Note 
that /C is a priori a generic vector of grade (—1), i.e. it lies in the same Lagrangian subspace 
as the vector Hq above, with the additional restriction of a vanishing component along R, 
since {R*,dR*) = by definition. However, we should note that the constant part of /C 
is not physical and can be absorbed into Rq by imposing the boundary condition that /C 
vanishes asymptotically. This choice is useful in the discussion of explicit solutions. 

We now turn to the Poisson equation ( [4. 441 ) for the vector T-Lq. As shown in ( 4.43| ), 



the source term of is along the unique grade (+3) component, that is along the vector R. 
We have therefore identified T-Lq as a vector lying in a Lagrangian submanifold containing 
R, whose components along the grade (—1) directions are n„ harmonic functions and the 



^^The non-generic solutions to this equation require constraints on the vector of T-duahty parameters 
and only exist if l4,('Ho,'Ho,'Ho,R) = 0, which would lead to solutions with irregular horizons. 
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component along the direction of i? is a single non-harmonic function, V. One can directly 
compute the source term by varying the combination in ( [4.29 ) and using ( 4.35| ), to find 



dT+ A *dT+no 



-(iT+ A*dT+[2Im(e-^" 






(4.47) 



which is explicitly proportional to the vector R, up to a real function. We can write this 
result in a simpler form by taking the inner product of ( 4.36| ) with R* and comparing with 



( I4.23D , to obtain 



v = -{R*,n) = -{R*o,no), 



(4.48) 



where we used ( |4.15| ) and (4.34). The Poisson equation for this function can now be found 
by projecting ( [4. 441 ) along Rq, as 



dirdV = -{R*Q, dT+ A i^dT+no) = -^ h{d)C, *d)C, Uq^R). 



(4.49) 



Note that, by ( 4.47 ), the right hand side of the first equality is nonvanishing, since it is the 
inner product of R with its magnetic dual. The second equality is a direct consequence of 



( 2.40 ). We moreover record the following relations for the expressions involving T-duality 
matrices in ( [4.29| ) and ( [4.47|) , 



dT+[2e"^Im(e" 
(iT+ A*dT+[2Im(e" 



V)- 


-f« 


V)- 


-^^ 



^exp(T+)/l(d/C,?^o,iZ), 
^hidlC,^dlC,no,R)R, 



(4.50) 



which are direct consequences of ( p. 351) and ( |2.40| ), in combination with ( [4.35 ). These 
equations, along with ( 4.22| ), allow us to evaluate the action of T-duality generators on the 
various objects relevant to the system. As the action of a finite T-duality is expressed as 
a finite sum due to the nilpotency of T+, it follows that one may also compute the action 
of finite T-dualities, as alluded to above. This is especially important when trying to find 



explicit solutions to this system [26]. 

The final quantities to be fixed are the angular momentum vector, oj, and the function, 
M, appearing in the expression for the scalars. In the single centre class, these are both 
harmonic and are dual to each other. Indeed, in that case the T-duality parameters in ( 4.1^ ) 
must vanish and the first of these equations imposes exactly that M and uj are harmonic. 
In the more general multi-centre case, one has to compute the nontrivial 7+ to obtain the 
analogous equation. This can be done straightforwardly using ( 4.18 ), ( 4.34D -( ^.36| ) and 



( 4.50 ) to show that 

-kdoj -dM = (rC^dw) 



{-Ho, dno - dT+ no) = 4 hidJC, ?^o, ^0, R) • 



(4.51) 



One can alternatively obtain the same expression by manipulating the definition of the 
non-harmonic function M in ([4.4|) using the above results. The integrability condition on 



(4.51) leads to the following Poisson equation for M: 

d*dM= -iCHo, dT+ * d-Ho) = ^ h{dlC, *dno,'Ho, R) , 



(4.52) 
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which can be solved once the grade —1 component of Tio and /C are chosen. 

This concludes our duality covariant presentation of the composite non-BPS system in 
terms of the real basis. In the next section, we summarise the final form of the equations 
to be solved and we comment on some properties of the solutions. 

4.5 Summary of results 

In this short section, we summarise all relevant formulae for the composite non-BPS system 
in the real basis. All relations presented here were shown explicitly in the previous sections 
and we refer to the discussion there for the details. We find it however useful for future 
applications to give a self-contained account of the final form of the system. 



The ansatze for the metric and gauge fields are given in ( |3.1[ ) and (3.3) in terms of 
the function e , the one- form lo and the spatial vector fields dw, while the electromagnetic 
potentials are fixed by ([4.5|). The first order flow equation for the composite non-BPS 
system is given by ( 4.30D , as 



i.dw = -[d-2 dT+] [2 Im(e-^-^° V) -^VR-fR*]. (4.53) 

Here, M, V, are functions to be specified below, while R and R* are a constant and a 
non-constant very small vector respectively, where {R,R*) = 4. The non-constant R* is 
related to a constant very small vector, Rq, by ( [4. 151 ) 

R* = exp[T+]i?S , (4.54) 

which also satisfies {R, Rq) = 4. In this and all equations in this section, Tj is a generator 
of the T-dualities leaving R invariant, parametrised by a vector of harmonic functions, /C. 
As discussed in section |2|, the vector of parameters /C lies in the grade (—1) component of 
the vector space according to the decomposition implied by the T-duality. It is therefore a 
three-charge vector satisfying 

h'^iR,R^,IC) = -{R,R*o)}C, (4.55) 

which indeed specifies a vector of n^ degrees of freedom. 

The solutions to the flow equation ( 4.53| ) are simplified by introducing a vector, Hq, 
of grade (—1) (+3), i.e. satisfying 

hi{R,R*Q,no) = -{R,Ro)no + 3{no,Ro)R. (4.56) 

Note that ( |4.55p is trivially a solution of the first equation, found by setting the grade (+3) 
component, (i?g,^o)) to vanish. The equations resulting from ( [4.53| ) take the form ( 4.36| ), 



2 e-^Im(e-^"V) = - exp[T+] (no - ^V R - f R*oj , (4.57) 

i.dw = exp[T+] {dUo - dTpio) , (4.58) 
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where V is now identified with the grade (+3) component of Ho, as F = {HqtRq)- The 



compatibihty relation for the last relations leads to the field equation for T-Lq, as in (4.44) 
and ( pO| ) 



dird-Ho = dT+Ai.dT+no = -^ hidJC^i.dlC^'Ho, R) R . (4.59) 

As the right hand side of this relation is only along R, it follows that all grade (—1) 
components of T-Lq are harmonic, whereas V is not, leading to ( [4.49| ), as 

di^dV = -^h{dlC,i<dlC,no,R), (4.60) 

by taking the inner product of ( |4.59| ) with Rq. The final dynamical equation required is 
the one for the function M in ( ^.58| ) and the angular momentum vector to, both of which 
are conveniently given by ( 4.5l| ), as 



•dw -dM = {no, dUo - dT+ no) = ^ hidlC, ?^o, ^o, R) ■ (4.61) 



Taking the divergence of this equation, one obtains a Poisson equation for M. 

These equations can be seen to be equivalent to the known formulation of the composite 
non-BPS system, as given in a fixed duality frame in [11, ^], by making a choice for the 



constant vectors R and Rq. In fact, these two papers use different frames for describing 
the system, both of which can be reached from our general formulation. The frame of |12| 
is found by choosing 

i?oc(0,5?), R*o^{6i,0), (4.62) 

where we disregard the (arbitrary) normalisation. Similarly the frame in [^] is found 



by interchanging the expressions for the two vectors in (4.62), corresponding to an elec- 
tric/magnetic duality. 

We close with some comments on the structure of the solutions. First, the physical 
scalars and the metric scale factor can be obtained by solving ( 4.57| ) in the standard way 



Since all quantities above are appropriate combinations of the single centre solution 



in 1 13], up to overall T-dualities, it is possible to use many of the results given there. For 



instance, the metric scale factor is given by 

e-''' = h{no-^VR-fR*o), (4.63) 

where we used the fact that the quartic invariant is by definition invariant under all T- 
dualities. This expression is identical to the corresponding one for the single centre class, 
despite the fact that the functions M and V are not harmonic in the present context. As 
in ||l^, one can simplify the expression for e as follows. The decomposition of the vector 
^0 — ^^ R ™ grades (—1) ® (+3) implies that Ii{no — \y R) is linear in V, as any other 
power of the grade (+3) would vanish identically. In particular we find that 

no-lvke^K^^f-'^ => h{no-lVR)=0, (4.64) 

so that we can show the equality Ii{no — \y R) = —hino)- It then follows that 

e-'^ =h{no -\VR)-^ {Rl,no - \VRf 

= -h{no)-M^. (4.65) 
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This equation implies that the ny grade (—1) harmonic components in T-Lq must correspond 
to a rank three charge (i.e. a large electric charge in five dimensions), so that ( [4.65| ) leads 
to a non-degenerate metric. 



Based on (4.65), we conclude that the system above indeed describes the interactions 
of black holes that are non-supersymmetric in isolation, since liil-io) must be negative 
globally and in particular at each centre, for a regular geometry. We refrain from giving an 
explicit expression for the physical scalars, as these would involve the action of an arbitrary 
abelian isometry with parameters /C on the physical scalars for the single centre class as 
given explicitly in |jl^. Of course, the scalars can be computed from the standard formulae 
in Q for any desired solution. 

The characteristic features of all solutions in this class is that each one of the centres 
must be of the non-BPS type in isolation, as explained above, and that the charges at all 
centres must commute with the vector R and must not commute with the vector Rq, by 
regularity. It is then clear that such solutions do not exist for generic non-BPS charges 
at the centres. In addition, once an allowed charge configuration is fixed, one cannot 
have arbitrary values for the moduli at infinity. This is due to the constraint on T-Lq^ 
which contains only n^ + 1 asymptotic constants and the fact that some of these constants, 
together with the 2ny parameters in R and i?g, turn out to parametrise charges that are 



not described by the poles of T-Lq, through ( 4.53 ) and ( 4.58 ). A simple example of this 



situation is given by the class of two centre solutions for models with n^ > 3. In this 
case, the two non-BPS charges are in fact arbitrary, since one can always find a choice 
of R and Rq for any pair of non-BPS charges. However, it is possible to show [^] that 
the asymptotic moduli are constrained to lie on a {n^ + 2)-dimensional hypersurface of the 
2n^-dimensional moduli space. Note that this is very different from the multi-centre BPS 
solutions, where solutions exist a priori everywhere in moduli space and walls of marginal 
stability arise only when the constrains implied by global regularity are imposed. In the 
present case however, local solutions seem to exist only in certain hypersurfaces of moduli 
space, while walls of marginal stability might still arise on these constrained surfaces. 

Finally, it is worthwhile commenting on the evaluation of T-dualities appearing as 
matrices in the equations above. As we have shown in several examples, one can avoid 
introducing explicit matrices, instead computing the action of T-dualities on the relevant 
vectors by use of the definition ( [4.46| ) and equations (|4.50| ). Indeed, we find that it is pos- 
sible to reduce all required computations to a recursive application of these three relations. 
As the grading involves only four subspaces, this procedure terminates after at most three 
steps. 



5. Almost-BPS system 

In this section we present in detail the characterisation of solutions to the almost-BPS sys- 
tem of equations, in analogy with the steps taken in the previous section for the composite 
non-BPS system. While the discussion here is self-contained, we will occasionally refer to 
section 0, in order to highlight similarities and recycle some results. The starting point is 
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the solution to the nilpotency condition given in ( |3.3£ ), which we repeat exphcitly here 



dU + -e'^^*duj 



1 



-e^df = e'^'Z'' + -n {Z + Af[n]n Zb - e'^'MMZ - e'^TlbZ^) , (5.1) 

where Z = Z{-kF), Za = Za{-*^F) are the central charges of the spatial field strengths F in 
( |3.3| ). The vector Qa is connected to a constant very small vector, R* , through ( p.41 )-( 3.43| ). 
For later reference, we also give the inverse relations, which read 



.e- ( dU + '-e^^ i.duj]+^ e'fMin] ^ ^ 



where we used as shorthand the one-form 
fi = Im 



e^M[n] ^dU + ^e'^^i.duj]+ e^Af[n] t f}"e„ id? 



(5.2) 



(5.3) 



The final required equation is the compatibility equation ( |3.40|) , which can be rearranged 



using (5^) and ( p7^ ) to obtain 
d{a - arg[Afp]]) = - 



^Im 



Y { dU + -e'^^ *duj 



-'"ff'efdf 



(5.4) 



Here, we used the function Y defined in the first equality of (2.50) and the fact that it has 
unit real part. Similarly, one can show using the flow equations above, that /x is also given 
by the expressions 



Re(e^M[ny- 



/^ 



e2^*da; + 2a[e-*"Z]) , 



: i e^^ duj - Im(e"^"a df) + i lYr^dY 



(5.5) 



that will be used in due time. 

In order to solve these equations, we follow a path similar to the last section, by 



considering the electromagnetic potentials in section |5.l| and using them to simplify the 
equations. The connection with T-dualities is shown in section ^.2| , while sections ^.3| and 
5.4 are devoted to the linear system of equations governing this system and its integration 
in terms of local functions respectively. The reader interested in applications can find a 



summary of the final form of the system in section 4.5 

5.1 The electromagnetic potentials 

As a first step towards the solution of the system, we decompose the field strength F in 
the electromagnetic potentials and the vector fields that define the conserved charges. The 
electromagnetic potentials for this system are computed by their definition (^) as 



dC = 2e2^K[Z(*F)V + Z*(*F)AV] 
= -2d {e^^e-'"V]) + ie^ fi R* 



(5.6) 
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where we used (|5.4| ) and the definition ( |3.41| ) to rearrange terms. The first term in (5^) is 



already a total derivative, so that the last term must combine into the derivative a vector, 
which is necessarily proportional to R* . Using the relation of R* to a constant vector in 
(3.43), this requirement leads to the condition 



e^ H = -2WRe[Z{R*)Zb{R*)e\df] + 2dW , (5.7) 

where we introduced the a priori arbitrary real function W . The result for the electromag- 
netic potentials takes the form 

C = -2 e^K[e-^"V] + WR\ (5.8) 

while the corresponding central charges are 

Z(C) = -ie^ e''' + W Z{R*) = -ie''^{e^ + iW e"^ Np]^) 
ZaiO = WZaiR*) = We'fMM^Qa. (5.9) 

We can now construct the vector potentials by the definition 

Z{dw) = Z{F) - Z{C) duj , Za{dw) = Za{F) - ZaiO du , (5.10) 

which leads to the expression 



•k dw = 2 Im 



d{e~^e-'''V) + {^e-^fi + Wdoj) 7^* + ie~>e~^AA[0]2 V ,(5.11 



where we use the shorthand in (15.31) and 

3ia ,^ ^ 1 



n* = e-— AA[0]2 V- e-^N[^]^n''DaV . (5.12) 



Here, it is worth pointing out that, unlike in the composite non-BPS system (cf. ( |4.9| )), 
there is no component of the vector fields that is vanishing a priori. Nevertheless, the 
projection of the vector fields in ( ^.11| ) along the available constant vector R* is still relevant 
and can be computed as 



Y\ \ \Y\ 



,-u\zm\J-uW ^^^^3, 



\Y\ \ \Y 

where Y is again as in ( |2.50| ). This can be simplified by imposing consistency of the two 
expressions for ^ in ( |5.5D with ( |5.7D , combined with ( |3.43D , leading to 

-uW\ Jrj\Z{W)\\ _^^_u\Z{W)\J^_uW 



l + 2e-^Jd e-^1^ =2e-^lj^d^e-^). (5.14) 



Now, it is simple to show that ( p.l3| ) and ( p. 14 ) imply that the projection {R* ,dw) is given 
by a harmonic function, V, defined as 

{R*,dw)= -2i.dle-^^-^^^]=-*dV, (5.15) 
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while the function W is fixed as 

l + 2e-^^ = /3y, (5.16) 

where /? is an arbitrary constant. Using these results, ( ^.7| ) simplifies to 

t^=\Y\^, (5.17) 

which will be used in due time. 

5.2 Connection to T-dualities 

We now discuss the relevance of T-dualities for the almost-BPS system, which will be 
important for the integration of the flow equations, as for the composite non-BPS system 
in the previous section. In order to exhibit this, we consider the very small vector 

R = Im [-M[n] V + n"" DaV] , (5.18) 

which is always mutually nonlocal with the constant vector R* and is not constant in 
general. As is clear by their definitions these two vectors are related in exactly the same 



way as the pair of very small vectors, R and R* in ( 2.17 ) and (|2.49| ) respectively, up to 
rescalings. Note that the situation is opposite to the one for the composite non-BPS system, 
where i? is a constant vector up to rescaling, while R* is not constant. 

Based on the discussion in section |2[ we associate R to the grade (-1-3) component of the 
decomposition ( p.29|) , while R* is identified as the corresponding grade (—3) component. 
Moreover, one can check that the normalised very small vector 

R = 2\Z{R*)\-^\Yflm[-M[n]V + n''DaV] , (5.19) 

has a constant inner product with R*, namely {R,R*) = 4. As explained in (^.5|)-(P^), 
the condition that i? is a very small vector can be generally written as 

R = exp[T-]Ro , (5.20) 

where Rq is a constant very small vector and exp[T~] is a T-duality matrix leaving R* 
invariant, parametrised by a grade (-1-1) vector of functions. It follows that the derivative 
of R can be expressed as 

dR-dT-R = 0, (5.21) 

which is closed by the property that T-dualities are abelian, exactly as in the composite 
non-BPS case. This is consistent with the known formulation of the almost-BPS system 



in five dimensions |14, 15, 16|, which can be written in terms of a T-duality parametrised 
by harmonic functions, acting on the scalar and vector fields. 

We can find the relevant T-duality parameters in ( ^.20D by explicitly computing the 
derivative of R, using the fiow equations ( ^.1| ) for the almost-BPS system. After a lengthy 
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but straightforward computation, we obtain that the derivative of ( 5.19 ) is indeed given 

by 



ZidR) =T-qZ{R) = Z{dT-R) , 
ZaidR) =T-QZa{R) = Za{dT-R) , 



(5.22) 



where we indicated that the result is given by the variation of R under the T-duahty 
transformation in ( p.62 ), as shown in ( |2.63D . The values for the one- form generators are 
given by 



7- 

Q-a 



\ [dU 



^Y e^^ * duj 



\ '^{Ye-^'^Vt ■ dt)] 



Y [cabc^^e^dt^ - e'^'eajdP + | J^a [e'^Q^ebjdt^ + 2 w e~'"njdt^ 



:^)]. (5.23) 



The expression for 7_ can be rewritten using ( |5.5[ ) in a form similar to the corresponding 
T-duality parameter for the composite non-BPS system in ( 4.18| ), as 



1 2U 
6 ^ 



{■kduj - dM) , 



(5.24) 



where we used the definition ( 2.50| ) for the function M in terms of the phases e*" and A/'[il]. 
In this form, it is manifest that setting the T-duality parameters j_ , Q°i to zero, one finds 
that the angular momentum is given in terms of a harmonic function, M, while the second 
of ( ^.231) becomes a reality constraint on the scalar flow, similar to ( 2.59 ) in the composite 
non-BPS system. This restriction therefore leads to the single centre subclass, which is 
common to both the composite non-BPS and almost-BPS systems. 

Applying the considerations of section |2| on the system at hand, we recall that a 
generic T-duality leaving R* invariant is parametrised by a rank three grade (—1) vector 
of parameters, which we denote by K, so that (|5.20| )-( ^l21 ) become 



R = exp[T^]i?o , dTj^Ro = dJC . 



(5.25) 



In what follows, we will generally suppress the explicit subscript /C from the T-duality 
generators for simplicity. One checks that /C is indeed a rank three vector, i.e. 



d.T'd.T^R 



h{dR) = hidT-R) = h{d)C) = , 
^ h{d^T-R, d,T-R, R*) = -±U{d^lC, d,IC, R*) 



(5.26) 

(5.27) 



These equations are clearly dual to the corresponding equations ( 4.22| ) for the composite 
non-BPS system. Alternatively, the same property follows from the fact that dT^Ro is of 
grade (-1-1) and the rank of such a vector is at most three. Indeed, one can directly verify 
that the vector dJC, as defined in ( ^.25|) , satisfies the constraint ( 2.64| ) and (|2.65| ), which 
are equivalent to the real constraint 

hi{Ro,R*,}C) = {Ro,R*))C. (5.28) 

These equations are explicit realisations of the general situation discussed in section 0. 
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Finally, we stress the difference between the constant very small vector Rq of grade 
(+3) and the original vector in ( 3.28 ), which was also used to define R in ( 5.1S| ). This 
can seem confusing, especially in view of the fact that we used the constant vector Rq to 
define the constraint in ( ^.28 ) and ultimately the grading of the vector space. However, 
as already explained in the analogous situation for the composite non-BPS system, below 
( |4.22 ), the grading associated to the T-dualities leaving invariant R* is only defined up to 
the action of the T-dualities themselves. In this respect, one can chose any constant vector 
-Ro in this orbit to define the grading. Given this redundancy in the definition, it will be 
convenient to fix T^ to vanish in the asymptotic region, such that Rq = R\r^oo- 

5.3 The linear system 

We are now in a position to use the above results to write the system of flow equations in 
the real basis, in terms of the symplectic section, V, the two very small constant vectors 
R* , Rq, and the relevant T-duality generators. 

To show this, we insert ( ^.16 ) and (^.17 ) in ( ^.11 ) to eliminate the spurious quantities 
W and fj, in favor of V and dV respectively. Moreover, it is useful to note the relation 



dT- [2 e-^ Im(e-'° V)] = 2 e~^ Im [-3 i 7* e"*" V + i 7* ^' A V + QI A V] , 



(5.29) 



as well as the identity 

{d-2dT-)VR = dVR-VdR, (5.30) 

which follows from ( 5.21| ), and we remind the reader that the explicit expression for dR is 
given by ( 5.22| ) with parameters as in ( 5.23| ). One can then verify that the expression 



i^dw = -{d-2dT- 



2Im(e 



-U—ia 



V) 



IVR 



M £)* 
■yR 



13 i.dojR\ 



(5.31) 



is equivalent to ( 5.11 ) above. This result is manifestly duality covariant in the real basis, 
as it is written in terms of real symplectic vectors only. In particular, note that ( ^.31| ) is 



completely analogous to the corresponding result ( 4.30 ) for the composite non-BPS system, 
up to the term explicitly proportional to the angular momentum. However, it is simple 
to shown that this term is unphysical, after considering the electromagnetic potentials as 
well. These are computed by combining the result ( |5.8D with ( 5.16D for the function W ^ to 
find 

C = -2e^5ft[e-*'^V] + (/3-^)i?*, 



which leads to the following expression for the total spatial field strengths 
F = Qduj + dw= -(2 e^K[e-*"V] + ^ R*\ du 
-^d-2dT') [2Im(e-^-*" 



(5.32) 



^V) 



VR 



M 
V 



R* 



(5.33) 



Given that the constant /3 does not appear in the gauge invariant total field strengths, we 
conclude it corresponds to a residual gauge transformation of the type A -^ A + 13 R* dt 
and therefore is unphysical. Henceforth we set /3 = in all relations, for simplicity. With 
this choice, ( |5.31 ) is formally exactly the same as its counterpart in the composite non-BPS 
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system in ( ^4.301 ) , up to changing the relevant T-dualities from those leaving R invariant to 
those leaving R* invariant. 

Note however that this choice of j3 is not the most convenient one for all purposes, 
as for example in showing that the above equations describe multi-centre BPS solutions 
as a particular case. It can be shown that this is the case when V is a constant, but one 
only recovers the standard form of BPS solutions after imposing /? = y, as is clear from 



equation (|5.32| ) 



5.4 Integration and local structure 

The presence of the T-duality connection dT" in ( |5.31| ) does not allow for a straightforward 
solution in terms of harmonic functions, but one can follow steps similar to the composite 
non-BPS system in order to solve the system in terms of local functions. We can write the 
scalar and vector fields as 

2e-^\iTi{e-''^V)-\VR-^R* = -H, dw = i^{d-2dT-)'H , (5.34) 

where the vector of functions % is the solution to the non-harmonic equation 

dirdn-2di.dT-'H-2dT- Ai.dn = 0. (5.35) 

This can be simplified and cast as a Poisson equation after introducing a rescaled vector, 
as 

Ho = expl-T-jn , (5.36) 

which in turn is the solution to the equation 

d • dHo - d * dT-Tio - dT" A •dT-'Ho = . (5.37) 

This can be formally obtained from ( [4. 371) upon exchange of T-duality transformations. In 
terms of the new vector, Tio, the scalar and vector fields are given by 



2e-^Im(e-*"V) = -exp[T-] t^Ho - ^V Rq - f R* 

dw = exp[T-] * dTio - *d(exp[T-]) no , (5.38) 

which is the final form of the system in the real basis. 

The solutions to the above system can be characterised using the fact that the compo- 



nents of the vector T-Lq are restricted, in the following way. The form (5.38) of the vector 
fields allows us to compute 

Za(exp[T-] * d-Ho) - Ar[n]na ^(exp[T-] • dno) 

= Za{*dw + exp[T-] dT-Ho) - J\f[n]na Z{i<dw + exp[T-] dT^T^o) 

= 2Q_a- e^'^Na + 2% e-^fta 7- (3 ^ - 2) . (5.39) 
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The crucial observation is that the central charges in ( 5.3S| ) satisfy the reality constraint 



( 12:591) , as 

e^" [Z^iexplT-jdno) - Af[n]n''Z{exp[T-]dno)] = 

c"*'=rJbZc(exp[T-] dUo) + n"" f Z(exp[T-] dUo) - J\f[n]n^ Zb{exp[T~] dUo) ) , (5.40) 



which restricts the components of the corresponding vector to lie on a particular Lagrangian 
subspace [|l3|. The same constraint holds for the integrated vector T-Lq, for which 



Za{eM^-]no) - Mima Z(exp[T-]?^o) 



o,e-^y. 



U^ 



(5.41) 



The real form of this relation is the same as for the composite non-BPS system ( 4.56| ), 
which we recall in this section for completeness 



-I^{Ro,R*,7io) — —{Ro,R*)T-Lo + 3{T-Lo,R*)Ro , {Ro,'Ho) — 



(5.42) 



where we used ( 5.25| ) to undo an overall T-duality on all terms in this equation. The 
constraints ( |5.42 ) are exactly dual to ( p.28p , as they are related by replacing R* and 
Rq. We therefore conclude that Hq and /C lie in opposite Lagrangian subspaces, i.e. the 
two vectors have no common directions and span 2n^ + 1 independent components in the 
2(nt, + l)-dimensional vector space. 

To be more precise, the vector /C is of grade (+1), as explained in section |5.2| , while 
the vector Tio and its derivative lie in the Lagrangian subspace composed by grade (—1) 



and (+3) components in the decomposition (2.29), as 

no G (E"")<-''©E<'' , 



(5.43) 



exactly as in ( 4.42| ). Applying this to (5.37), the following pattern arises for the various 
terms 



d • dTio -d* dT-Tio - dT' A ^dT'T^o = , 

(-i)ffi(+3) (-3)e(+i) (-1) 



(5.44) 



where we used the fact that T" lowers the grade of a vector by (—2). In direct correspon- 
dence with (f4.43| )-( 4.45 ) for the composite non-BPS system, we find that ( 5.44 ) decomposes 
into two equations according to its graded decomposition, as 



d • dTio - dT- A i^dT-Ho = ; 
d * dT- = , 



(5.45) 
(5.46) 



where in the second equation we used the property that no T-duality T leaves the vector 



( 5.43 ) invariant. This is true because the grade 3 component of T-Lq can be identified as the 
nowhere vanishing harmonic function V defined in (|5.15|), as 



{R*,no) = -v, 



(5.47) 
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as can be seen by contracting ( 5.38| ) by R*. We now analyse each of the two equations 
(|5l5D-(|5l^) in turn. 

The solution to ( |5.46| ) is equivalent to the condition 

dirdlC = 0, (5.48) 

where we used ( 5.25| ). It follows that the vector of parameters /C is a generic grade (+1) 
vector of harmonic functions, /C, which is of rank three. Note that, in this system, the poles 
of /C represent new independent physical charges, since this vector is by definition linearly 
independent from Hq. In fact, the poles of this function are the only relevant information. 



since one may always absorb the constant part of /C into Rq in (5.25). 

We now turn to the Poisson equation ( |5.45|) and observe that the source term is of 
grade (—1), according to ( ^.44 ). It follows that only Uy — 1 out of the n^ components of 
T-Lq are sourced, leading to an equal number of non-harmonic functions. The remaining 



component is the harmonic function V, already identified in ( |5.47 ) above. The source term 
in ( 5.45 ) can be computed explicitly using ( p.62| ), with the result 

dh 



Z (- exp[T-] dT- A i^dT-Uo]) 
Za (- exp[T-] dT" A i<dT-no) 



V^(R*,dR,*dR) 
64 dZ^ ' 

-.±V~{R*,dR,i<dR). 

64 az«^ ' 



(5.49) 



i.e. it is proportional to the vector defined in ( 5.27] ). This vector is of grade (—1) by 



construction and it can be verified to satisfy the constraint in (|5.40D . We can now rewrite 

(5.50) 



(15.451) as 

di.d'HQ = ^Vl'^{R*,dJC,i^dlC), 

where we used the definition of /C in ( ^.25 ) and ( |5.27] ). 

Finally, we present the covariant form for the equation determining the angular mo- 
mentum and the function M in (2.50) and ( |5.24| ). The starting point is the first of ( 5.23| ), 
which upon use of (|5.5| ) can be written as 

*doj - dM = {exp[T-]no,*dw) = CHq, dTio - dT" Hq) , (5.51) 

where we also used ( ^.38 ) in the second equality. Explicit computation of the last expression 
using the definition ( |2.36 ) leads to the alternative form 



i^duj — dM 



lV{no,dR) = -^V{'Ho,d)C). 



(5.52) 



2 ' \".u,"-v — 2 

Taking the divergence and the curl of this equation one obtains the relevant equations for 
the function M and the angular momentum respectively. The resulting Poisson equation 
for M reads 

di^dM = ld{V{no,i^dlC)), (5.53) 

and can be solved once Ho and /C are specified. Note that upon setting the parameters, /C, 
of the T-dualities to vanish, these equations imply that M is a harmonic function, while cj 
is the corresponding dual one-form, consistent with the single centre class. 

This concludes our duality covariant presentation of the almost-BPS system in terms 
of the real basis. In the next section, we summarise the final form of the equations to be 
solved and we comment on some of the properties of solutions. 
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5.5 Summary of results 

In this short section, we summarise the relevant formulae for the almost-BPS system in 
the real basis. All relations presented here were shown explicitly in the previous sections 
and we refer to the discussion there for further details. We find it however useful, both 
for clarity and for future applications, to give a as self-contained as possible account of the 
final form of the system. 

The ansatze for the metric and gauge fields are given in ( p. ID and ( p.3p in terms of 
the function e^, the one- form uj and the spatial vector fields dw, while the electromagnetic 
potentials are fixed by ( 5.32| ). The first order equation for the almost-BPS system is given 



by (lOl 



as 



i^dw = -[d-2 dT-^\ [2 Im(e-^ -*° V) -^VR-fR*]. (5.54) 

Here, M, V ^ are functions to be specified below, while R* and R are a constant and a 
non-constant very small vector respectively, where {R,R*) = 4. Here, the non-constant R 
is related to a constant very small vector, iJg, by ( ^.20| ) 

R = exp[T^] Ro , (5.55) 

which again satisfies {Rq,R*) = 4. In all equations, T^ is a generator of the T-dualities 
leaving R* invariant, parametrised by a vector of harmonic functions, /C. As discussed in 
section ^, the vector parameter /C lies in the grade (+1) component of the vector space 
according to the decomposition implied by the T-duality It is therefore a three-charge 
vector satisfying 

h'{Ro,R*,IC) = {Ro,R*)lC, (5.56) 

which indeed specifies a vector of riy degrees of freedom. 

The solutions to the flow equation ( 5.54| ) are simplifled by introducing a vector, I-Lq, 
of grade (—1) © (+3), i.e. satisfying 

hi{Ro,R*,no) = -{Ro,R*)no + 3{no,R*)Ro . (5.57) 



Note that ( ^.56 ) follows from a similar constraint, obtained by interchanging R* with Rq, 



that projects to the (+1) ® (—3) component of the vector space. The equations resulting 
from ( 5.54| ) upon use of Tio in ( |5.57| ), take the form ( ^.38| ), 



2 e-^Im(e-*"V) = - exp[T-] f ?^o - i ^ ^o - f R*) , (5.58) 

*dw = exp[T^] {drCo - dT^Tio) , (5.59) 

where V is now identifled with the grade (+3) component of Hq, as 1/ = {Ho,R*). The 
compatibility relation for these relations leads to the field equation for T-Lq, as in ( 5.50| ) 



d*d'Ho = dT^A*dT- = ^V l4{dJC,i.dlC,R*) . (5.60) 
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As the right hand side of this relation is a vector of grade (—1), the corresponding compo- 
nents of Tio are not harmonic, whereas V is, as can be seen by taking the inner product of 



( lOOD with R* 

{R*,di.dno) =di.dV = 0, (5.61) 

where we used (|5.56|). The final dynamical equation required is the one for the function 



M in ( 5.58| ) and the angular momentum vector oj, both of which are conveniently given by 



as 



•dw -dM = {Uq, dTio - dT^ 'Hq) = \v (/C, Uq) . (5.62) 

Taking the divergence of this equation, one obtains a Poisson equations for M. 

The equations above can be seen to be equivalent to the known formulation of the 
almost-BPS system, as given in five dimensional supergravity [|l^, 15, 16 1, by making 
choice for the constant vectors Rq and R* . Indeed, upon choosing 



a 



fiooc(0,5?), i?* a (5^,0), (5.63) 

where we disregard the (arbitrary) normalisation, one can show a complete equivalence 
of the above to the original system in |14|. This particular frame is convenient in that it 
allows to lift to five dimensional solutions that are locally but not globally supersymmet- 
ric. However, our formulation of the almost-BPS system is closed under four dimensional 



dualities and includes all duals of the system in |14]. More recently, it was shown in |29| 
that some of the EPS structure is preserved in four dimensions as well, upon reinterpreting 
the constant vector Rq as Fayet-Iliopoulos terms in a gauged theory. 

We close with some comments on the structure of the solutions. First, the physical 
scalars and the metric scale factor can be obtained by solving ( 5.58| ) in the standard way 



once T-Iq and M are solved for. Since all quantities above are appropriate combinations 



of the single centre solution in [13|, up to overall T-dualities, it is possible to use many of 



the results given there. For instance, the metric scale factor is given by (4.65), as 

g-4C/ ^ -i^(ji^^) - m2 , (5.64) 

in exactly the same way as for the composite non-BPS system. However, in this case the 
situation is richer and more complicated, in view of the fact that the grade (—1) components 
of T-Lq are not harmonic and its grade 3 component V does not necessarily carry a pole 
at all centres. Indeed, it turns out that not all black holes described by the almost-BPS 
system are non-supersymmetric in isolation. On the contrary, the presence of both BPS 
and non-BPS types of centres, is the distinguishing property of this system, as shown in 
1 14, 15, [l^. Clearly, the fact that the harmonic functions /C lie in a subspace independent 



of the one where 1-Lq lives is the crucial ingredient that allows for both BPS and non-BPS 
types of charges to exist simultaneously. 

As seen in the case of the composite non-BPS system, solutions do not exist for all 
charge configurations and this holds also in the almost-BPS system. Moreover, it is not 
possible to obtain arbitrary asymptotic moduli for a given allowed charge configuration, for 
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exactly the same reasons explained in section 4^. Indeed, ( [4.53 ) and ( [4.58| ) have exactly 



the same structure in both cases, so that some of the n^ + 1 asymptotic constants in T-Lq and 
the parameters of R and Rq will correspond to charges rather than moduli. We once again 
refer to [^ for more details on the structure of almost-BPS solutions in four dimensions 
and for explicit examples. 

6. Conclusion 

In this paper, we gave a comprehensive treatment of the flow equations describing multi- 
centre under-rotating black holes mAf = 2,D = 4 supergravity coupled to vector multiplets 
with a symmetric scalar manifold. In particular, we considered the non-linear sigma model 
obtained after timelike dimensional reduction to three dimensions and derived the general, 
frame independent, flow equations for two systems of multi-centre non-BPS black holes, 
namely the composite non-BPS and almost-BPS systems. 

This represents a generalisation of the systems given in specific frames in |[ll|, |lj], to 
systems that are closed under electric/magnetic duality. The resulting structure for the 
vector fields and scalars in terms of real symplectic vectors turns out to be very similar 
for both systems. In particular, both systems are described in terms of space-dependent 
transformations along abelian subgroups of isometrics on the scalar target space. In terms 
of the natural embedding to string theories, these subgroups of the full duality group 
are conjugate to the so-called spectral flow transformations, that are combinations of T- 
dualities with gauge transformations on the p-form gauge fields. In this paper, we refer to 
them simply as T-dualities for brevity. 

The main distinction between these solutions and the BPS multi-centre solutions, is 
that the electromagnetic vector fields are not harmonic anymore, but satisfy instead 

di^dw-2dT^ Ai^dw = , (6.1) 

where the functions /C are themselves harmonic. Note that the consistency of this equation 
requires that the generators T are indeed abelian, as for T-dualities. It follows from 
this equation that the poles of the harmonic functions /C contribute to the electromagnetic 
charges in a non-linear way. Despite the interpretation of these functions as parameters 
of abelian isometrics of the scalar manifold, they are not associated to a gauging of the 
theory. 

The crucial property that makes a general discussion in terms of covariant objects 
possible is that the action of general T-dualities can be given explicitly using the quartic 



invariant, I4, of symmetric special Kahler geometry. Indeed, as summarised in sections 4.5 



and 5.5, all relevant equations are written in terms of this invariant only, evaluated for the 
real vectors parametrising the solutions. In this form, these duality covariant systems are 
not significantly more complicated than the equations given for the composite non-BPS 



system in |11] and for the almost-BPS system in [14| and can be solved in exactly the same 
way. 

The main advantage of the formulation displayed in this paper is that one need not 
define solutions in a fixed duality frame in terms of generic parameters, and only compute 
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the electromagnetic charges and asymptotic moduh a posteriori, as in the constructions 
of |l^, ^ |ll|, |l^. In contrast, one can start from any configuration of physical charges 
satisfying the required criteria associated to each system and construct the corresponding 



solution, using the results summarised in sections 4.5 and 4.5. This is in particular very 
useful for studying the domain of stability of these solutions in moduli space. Using this 
formulation, one can start from a given set of electromagnetic charges consistent with the 
system (e.g. they have to all mutually commute with a common very small vector R in 
the composite non-BPS system), and parametrize the most general very small vectors R 
and R* satisfying the corresponding constraints. Using the formulation of this paper, one 
can then determine the most general solution associated to a given charge configuration 
and define the domain of existence of such solutions in moduli space. As opposed to BPS 
solutions, the domain of existence of such solutions in moduli space will be restricted to 
a hypersurface of non-zero co-dimension. The normal directions to the hypersurface are 
probably not forbidden physically, but rather push us out of the domain where we know how 
to describe the solution. For example, the BPS solutions within the almost-BPS system 
only exist on a co-dimension one hypersurface, but only due to the fact that the charges 
at all centers must be compatible with a single constant vector Rq, leading to a subset of 
all BPS solutions. Nonetheless, one can still wonder if there are walls of marginal stability 
within the hypersurfaces defined by each of the two systems described in this paper, i.e. 
whether there are boundaries of the domain of existence of such solutions at finite values 
of the moduli. We intend to study two-centre configurations in the aim of exhibiting (or 
not) walls of marginality for non-BPS solutions in a forthcoming publication. 

From a more general point of view, the unified description of the two known non-BPS 
systems and its relative simplicity are encouraging for further uncovering the structure 
of non-BPS solutions in supergravity. In particular, the isometrics of the scalar manifold 
seem to play a crucial role not only in the effective three-dimensional theory, but also in 
the real formulation in four dimensions. It would be interesting to understand the role of 
these isometrics in the reduction of the equations of motion to first order systems, which is 
not clear from our treatment in terms of nilpotent orbits. In fact, it is known that higher 
orbits, describing more complicated systems of non-BPS solutions, exist and one might 
hope that similar structures as the ones described in this paper appear in those cases as 
well. 
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A. A/^ = 2 supergravity and symmetric special Kahler geometry 

The bosonic Lagrangian of A/" = 2 supergravity coupled to n^ vector multiplets reads 
m, |6| 

Svr e-i C = -\R-i {D^V, D^V) - \ F^^ Gf . (A.l) 

Here, the -F^^, = d^Al, — d^Aj^ for / = 0, . . . n^ encompass the graviphoton and the gauge 
fields of the vector multiplets and G^"^ are the dual field strengths, defined in terms of 
the F^j^ though the scalar dependent couplings, whose explicit form will not be relevant in 
what follows. The gauge field equations of motion and Bianchi identities can then be cast 
as a Bianchi identity on the symplectic vector 

whose integral over any two-cycle defines the associated electromagnetic charges through 

(A.3) 




The physical scalar fields f, which parametrize a special Kahler space M4 of complex 
dimension TT-t, , only appear in ( |A.1| ) through the section, V, of aholomorphic [/(l)xS'p(2n^ + 
2, R) bundle over AI4. Choosing a basis, this section can be written in components in terms 
of scalars X^ as 

where F is a holomorphic function of degree two, called the prepotential, which we will 
always consider to be cubic 

F- K x^x^x^ _ M[X] 

for completely symmetric Cjjfc, i = 1, . . . n^, and we introduced the cubic norm A/'[A]. The 
section V is subject to the constraint 

(V,V)=i, (A.6) 

and is uniquely determined by the physical scalar fields t* = ^ up to a local f/(l) trans- 



formation. The U{\) gauge invariance of (A.l) is ensured by the appearance of the Kahler 
connection Q^j in the covariant derivative. The Kahler potential on AI4 is defined up to 
an arbitrary holomorphic function f{t) as 

K = -h^{iN[t-t\)+f{t) + f{t) (A.7) 

and we fixed the U{\) gauge invariance in terms of Kahler transformations by requiring 
that the Kahler connection is determined by the Kahler potential as 

Q = ^diJCdf] , (A.8) 
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such that 

gij = didj}C, D^V = {d^ + iQ^)V = D^Vd,,f = {diV + ^d^}CV)d^t\ (A.9) 

where DiV is the corresponding Kahler covariant derivative on the components of the 
section. With the prepotential ( A.5| ), the special geometry identities |37| reduce to 



which are used throughout the main text. 

We introduce the following notation for any symplectic vector J 



(A.IO) 



Z(J) = (J,V) 



Zi{J) = {J,DiV), 



(A.ll) 



with the understanding that when the argument is form valued, the operation is applied 
component wise. For instance, the central charge of the gauge field is 



ZiT) = ef (Go + fGi + ^CijktH^F'' - M[t]F' 



(A.12) 



for the prepotential ( A. 5 ). With these definitions it is possible to introduce a scalar de- 
pendent complex basis for symplectic vectors, given by (V, DiV), so that any vector J can 
be expanded as 

J = 29[-Z( J) V + g'WjZ{J) DjV] , (A.13) 

whereas the symplectic inner product can be expressed as 

(Jl, J2) = 29[-Z(Ji) Z(J2) + Za{Jl) Z%J2)] . (A.14) 



Finally, we introduce the notion of complex selfduality of the gauge fields ( A. 21 ), which 
satisfy the identity 

Jjr=_*jr^ (A.15) 

where J is a scalar dependent complex structure defined as 



JV 



-iV , JAV = iDiV . 



(A.16) 
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